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SECTION 1.0 


INTRODUCTION 


One of the more important recent developments in gas turbine blade materials 
has been the introduction of directionally solidified and single crystal 
castings. This casting process has matured to the level where it is now 
routinely used in the production of commercial and military aircraft jet 
engine turbine blades. Among the advantages of these materials are the 
following: 

o By eliminating grain boundaries, the high temperature creep strength 
is substantially increased. 

o The elastic modulus in the direction of solidification is lowered, 
reducing thermally induced stresses and, thereby, increasing low 
cycle fatigue life. 

o Corrosion and oxidation resistance are enhanced by the absence of 
grain boundaries. 

o Eliminating grain boundary strenthening elements raises the melting 
temperature and permits greater heat treatment flexibility. 

Unfortunately, these metallurgical and processing advances have not been 
matched by correspondi ng advancements in the knowledge and understanding of 
the mechanics of these materials, their failure mechanisms, and methods for 
life predictions. In order to realize the full potential of these materials, 
it is necessary to have a complete knowledge of the full envelope of life 
limiting parameters. Anisotropy introduces many life prediction questions 
especially for stresses which are not parallel to the direction of 
solidification. Oxidation resistant coating adds another layer of complexity 
to these questions. All of these issues are addressed in this HOST program. 

The program is divided into a base program which covers a duration of 
forty- two months and two optional programs which are to be exercised at the 
discretion of NASA. In the base program a primary and alternative coated 
single crystal material , operating at relevant turbine airfoil temperatures, 
are being investigated. In Option 1 the same two single crystal materials, in 
an uncoated condition and operating at root attachment temperatures and 
notched conditions, will be studied. In Option 2 a directionally solidified or 
recrystallized material, in a coated and uncoated condition, will be studied 
at temperatures occurring at the airfoil and root attachment, respectively. 

In the base and optional programs, candidate constitutive and life prediction 
models are being developed concurrent with laboratory specimens which will be 
tested under a variety of mechanical and thermal load histories to provide 
data for the final model selections. The models will be incorporated into a 
computer code which will be checked for operability on a representative 
turbine blade section. 
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The first year effort of the contract has involved material selections 
specimen preparation, basic material tests, literature searches of appropriate 
constitutive and life prediction models, initial formulation of constitutive 
models, and initial constitutive and fatigue tests. The specific 
accomplishments are listed in Section 5.0 and are expanded upon throughout 
this report. 

During the second year, the Task III, Level I, experiments will be completed 
on the primary single crystal alloy and the two coatings. Metallurgical 
evaluation of the results, together with constitutive and life prediction 
modeling, will continue in parallel with these tests. 
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SECTION 2.0 


TASK I - MATERIAL/COATING SELECTION AND ACQUISITION 


2.1 Material Selection 

A survey was conducted of the single crystal materials and surface protection 
coating systems currently in use or in an advanced state of development for 
use in turbine blade applications. Two single crystal materials and two 
coating systems have been selected for this program. 

2.1.1 Single Crystal Airfoil Material 

The trend toward improved engine efficiency and improved durability has placed 
increasing demands on materials for turbine airfoils that operate in the hot 
section of the gas turbine engine. To meet these demands, nickel -base 
superalloys have been designed specifically for use in single crystal form, 
which offers greater creep fatigue resistance at elevated temperatures and 
reduced susceptibility to corrosion or oxidation than conventionally cast 
alloys. Table 2-1 lists some of the single crystal alloys that have been 
developed for gas turbine engines. These alloys do not contain grain boundary 
strengthening elements such as carbon, boron, and zirconium; therefore, the 
incipient melting temperatures are higher than those of conventionally cast 
alloys. This allows for a higher temperature solution heat treatment, which, 
in turn, further improves alloy strength by providing a high volume fraction 
(approximately 65 percent) of fine, cuboidal, homogeneously distributed gamma 
prime precipitate. 


Table 2-1 

Single Crystal Superalloys Used in Gas Turbine Industry 
Alloy Composition (Weight Percent) 



la 

W 

Mo 

Cr 

A1 

Ti_ 

Cb 

Co 

Hf 

C 

B Zr Ni_ 

PWA 1480 

12 

4 

-- 

10 

5 

1.5 

__ 

5 

-- 

-- 

— — Bal* 

Rene N4 

4 

6 

1.5 

9.2 

3.7 

4.2 

0.5 

7.5 

— 

— 

- - Bal* 

CM SX-2 

6 

8 

0.6 

8 

5.6 

1 

__ 

4.6 

— 

-- 

— — Bal* 

CM SX-3 

6 

8 

0.6 

8 

5.6 

1 

— 

4.6 

0.1 

— 

- ~ Bal* 

*Balance 
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PWA 1480 was selected as the primary single crystal superalloy. It was the 
first superalloy designed specifically for use in single crystal form and was 
developed with the goal of achieving an optimum balance of creep strength, 
thermal fatigue strength, and oxidation and hot corrosion resistance. PWA 1480 
is the most widely used single crystal alloy in gas turbine engines today and 
the most advanced turbine airfoil material utilized in Pratt & Whitney 
production engines. PWA 1480 was certified for commercial use in the 
JT9D-7R4D/E engine in late 1981 and has since been certified for use in the 
JT9D-7R4G/H and PW2037 engines. 

Two heats of PWA 1480 have been procured for this program from the Howmet 
Turbine Components Corporation, Alloy Division, Dover, New Jersey. The primary 
heat, identified by Howmet as 200A14824, has been designated P9866. The 
secondary heat, identified by Howmet as 200B14773, has been designated P9867. 
The chemical compositions for the heats are compared to the nominal 
specifications in Table 2-II. The typical solution heat treated microstructure 
is illustrated in Figure 2-1 (A). 


Table 2-II 

Single Crystal Superalloys to be Studied in Contract 

Alloy Composition (Weight Percent) 


Alloy 

Heat 

Code 


Cr 

Co 

Elements 

11 £L 

li 

W 

Mo 

C 

PWA 1480 

Nomi nal 

Bal* 

10.0 

5.0 

1.5 

5.0 

12.0 

4.0 

— 

— 


P9866 
(Heat A) 

Bal* 

10.35 

5.5 

1.44 

4.95 

12.2 

3.9 

— 

0.01 


P9867 
(Heat B) 

Bal* 

10.3 

5.3 

1.44 

4.9 

10.3 

4.0 

— 

0.004 

Alloy 185 

Nomi nal 

Bal* 

— 

— 

— 

6.8 

— 

6.0 

14.0 

0.04 


P9921 

Bal* 

- «. 

-- 

0.001 

6.82 


6.10 

13.85 

0.04 


♦Balance 
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Typical Solution Heat Treated Microstructure Illustrating. 
Gamma /Gamma Prime Eutectic Islands in Gamma Matrix With Fine 
Unresolved Gamma Prime Precipitates of: A) PWA 1480, and B) 
Alloy 185. ( 500X Mag., Etchant: Mixed Acids) 
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Alloy 185 was selected as the secondary single crystal alloy with the intent 
of testing the range of applicability of the constitutive and life prediction 
models developed in this program. Alloy 185 has a higher volume fraction 
(approximately 70 percent) of gamma prime compared to PWA 1480 due to a higher 
aluminum content. In addition, the high molybdenum content and high volume 
fraction of gamma prime produce a morphology that is semi -continuous rather 
than cuboidal. As a result. Alloy 185 exhibits greater creep anisotropy than 
PWA 1480. 


A single heat of Alloy 185 has been procured for this program from the Howmet 
Turbine Components Corporation, Alloy Division, Dover, New Jersey. This heat, 
identified by Howmet as 242A15847, has been designated P9921 . The chemical 
composition for this heat is compared to the nominal specifications in Table 
2-1 I. The typical solution heat treated microstructure is illustrated in 
Figure 2-1 (B). Tensile and creep properties for this alloy are given in 
Figures 2-2 and 2-3. 




«3 
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CO 
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Figure 2-2 Tensile Properties of Cast Single Crystal Alloy 185 


2-4 


552 



Figure 2-3 500-Hour Rupture Life of Cast Single Crystal Alloy 185 


2.1.2 Coatings 

For gas turbine applications of high temperature alloys, the use of protective 
coatings is required to achieve practical lives for components which are 
designated to operate in oxidizing environments at elevated temperatures. One 
overlay and one diffusion coating were selected for this program to be 
representative of those employed in actual gas turbine engines. 

PWA 286 was selected as the overlay coating. It offers enhanced oxidation and 
corrosion resistance for turbine blades and vanes. The composition and 
deposition process are provided in Table 2-1 II . The typical microstructure for 
this coating is illustrated in Figure 2-4(A). 
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Table 2-1 II 


Coating Compositions and Processes 


Coati ng Type 

PWA 286 Overlay 


Composi tion 
NiCoCrAl Y+Si+Hf 


Deposi ti on 
Process 


Vacuum Plasma Spray 


PWA 273 Aluminide Ni A1 Pack Cementation 

(Outward 
Diffusion) 


PWA 273 was selected as the diffusion coating. It is an outward diffusion 
aluminide coating for nickel -base turbine blades and vanes. The composition 
and deposition process are given in Table 2-1 II. The typical microstructure 
for this coating is illustrated in Figure 2-4 ( B ) . 

2.2 Specimen Preparation 

2.2.1 Casting 

The superalloys employed for this program have been cast using the directional 
solidification process which is schematically illustrated in Figure 2-5. The 
molten metal, initially freezing on the chill plate, has a very fine 
polycrystalline structure. As the solidification front moves vertically, 
elongated grains (or columns) are formed. The dynamics of solidification are 
such that grains with a [001] crystal orientation parallel to the direction of 
heat flux grow preferentially compared to grains with other orientations. To 
produce a single crystal casting, a helix is used to act as a filter 
constricting the growth of the mar*y columnar grains that grow from the water 
cooled chill plate. The helix permits a single <001 > oriented grain to pass 
through. For orientations other than <001> (i.e., <111>, <11 0>, and <123>), a 
seed crystal of the desired orientation is employed in the mold to override 
the natural tendency for growth in the <001> direction. When molten alloy is 
poured into the mold, the seed crystal acts as a starting block for 
solidification. Directional solidification proceeds with the crystallographic 
growth direction dictated by the orientation of the seed crystal producing a 
single crystal casting of the desired orientation. 

Cylindrical single crystal bars with nominal 15.2 cm (6.0 in) lengths and 2.54 
cm (1.0 in) and 1.59 cm (0.625 in) diameters have been cast. The standard 
<100> crystallographic direction along with the <111>, <11 0>, and <123> 
directions have been chosen as the primary orientations for single crystal 
casti ng. 
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OF POOR QUALITY 



Figure 2-4 


Typical Micrographs of: A) PWA 286 Overlay Coating, and B) PWA 
273 Diffusion Coating Illustrating the Microstructural 
Differences Between the Coatings. Note the small interdiffusion 
zone associated with the overlay coating compared to that of the 
diffusion coating. The substrate is PWA 1480. (500X Mag., 
Etchant: Mixed Acids) 






Figure 2-5 Schematic Representation of Single Crystal Casting Process 


After casting, the bars were solution heat treated employing a 
step-temperature cycle with a thirty-minute hold at the 1293°C (2360°F) 
maximum temperature, followed by a forced gas cool to obtain the desired 
microstructure of fine gamma prime precipitate homogeneously distributed in 
the gamma matrix without the onset of incipient melting. The following 
inspection procedure was employed to ensure the quality of the castings prior 
to specimen fabrication: 

1) Grain etch - visual inspection for surface defects. 

2) Laue back reflection x-ray diffraction to determine the deviation 
from the primary crystal! ographic orientation and the presence of any 
low angle grain boundaries. 

3) Metal 1 ographic examination of the heat treated microstructure to 
ensure complete gamma prime solutioning without incipient melting or 
excessive porosity. 

Table 2-IV lists the matrix for the number of acceptable bars cast to date. 
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Table 2-1 V 


Matrix for the Number of Acceptable Cast Bars Based on 
Crystallographic Orientation and Bar Diameter 


Number of 2.54 cm (1.0 in) Number of 1.59 cm (0.625 in) 

Diameter Bars Diameter Bars 


Alloy 

<100> 

<111> 

<110> 

<123> 

<100> 

<111> 

<11 0> 

<123> 

PWA 1480 

82 

28 

73 

69 

53 

42 

37 

31 

(P9866) 
PWA 1480 

28 

3 

3 



8 



(P9867) 
Alloy 185 

15 




15 





2.2.2 Specimen Fabrication 

Specimen geometries were chosen to allow conditions comparable to those found 
in actual turbine blades to be produced during testing. PWA 1480 specimens 
(and, in the future, Alloy 185 specimens) were machined from 2.54 cm (1.0 in) 
and 1.59 cm (0.625 in) diameter bars. For coated property testing, the 
specimens were coated using the standard methods as listed in Table 2-1 I I . 
Specimens for determining coating material properties of the overlay coating 
were fabricated from ingots of hot isostatically pressed (HIP) powder. 
Diffusion aluminide coating specimens were prepared by applying the coating to 
specially designed specimens of PWA 1480. Additional details on the specimens 
are as follows: 

o Physical Property Specimens - The determination of the thermal 

conductivity, thermal expansion, and specific heat for the materials 
chosen for this project has been contracted to the Southern Research 
Institute, Birmingham, Alabama. Figure 2-6 illustrates the specimen 
geometries required for two of the physical property tests. The 
specification for the specific heat specimen is that a minimum volume 
of 16.4 cubic cm (1.0 cubic in) of material be supplied. Six 
specimens (two for each property) of each material are to be tested. 

o Single Crystal Alloy and Single Crystal Alloy /Casting Mechanical 

Property Specimens - Figure 2-7 (A) and 2-7 ( B ) illustrate the specimen 
geometries employed for coated and uncoated tensile and creep 
testing, and uncoated cyclic constitutive testing. Figure 2-7(0 is a 
diagram of the current standard strain-controlled thermomechanical 
fatigue (TMF) specimen used at Pratt & Whitney. The specimen was 
employed to determine the feasibility of using external extensometry 
(which would eliminate the need for internal ridges in the specimen) 
instead of the standard internal extensometry by performing direct 
comparison TMF testing trials. Table 2-V lists the number of 
specimens prepared to date. 
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A) THERMAL CONDUCTIVITY SPECIMEN 



B) THERMAL EXPANSION SPECIMEN 



Figure 2-6 Specimen Designs for Determining Physical Properties 


DIA. 
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A) TENSILE AND CREEP TEST SPECIMEN 


Figure 2-7 



B) CYCLIC CONSTITUTIVE TEST SPECIMEN 



C) INITIAL THERMOMECHANICAL FATIGUE SPECIMEN 



Specimen Designs for Single Crystal Alloy and Single Crystal 
Alloy /Coating Mechanical Property Tests 
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Table 2-V 


Number of PWA 1480 Mechanical Property Test 
Specimens Fabricated to Date 





Number of 

Specimens 


Type of Specimen 

Coating 

<100> 

<m> 

<110> <123> 

Tensile and Creep 
Test Specimens 

Uncoated 

14 

12 

14 

14 

PWA 286 

10 

9 

— 

— 


PWA 273 

10 

9 

— 

— 

Cyclic Constitutive 

Uncoated 

11 

11 

11 

11 

Test Specimens 






Thermomechanical Fatigue 

PWA 286 

12 

a. 

■»« 

cmmm 

Test Specimens 

PWA 273 

12 





o Coated Material Mechanical Property Specimens- Figure 2-8 illustrates 
the specimen geometries employed for testing the mechanical 
properties of bulk PWA 286 overlay coating material. The specimen 
diagrammed in Figure 2-8(A) was machined from PWA 286 ingots of hot 
isostatically pressed (HIP) powder. Figure 2-8(B) illustrates 
specimens fabricated from thick sheets of plasma sprayed PWA 286. The 
thick sheets were produced of plasma spraying thick layers of PWA 286 
onto substrates. The substrates were subsequently removed by 
machining. 

Figure 2-9 illustrates the specimens employed for determining coating 
material properties when applied to PWA 1480 substrates. Figure 
2—9 ( A ) diagrams the substrate for the PWA 273 coated specimen for 
tensile and creep testing. Two coating thicknesses will be used. The 
reasons for this particular configuration selection is discussed in 
Section 4.1. Figure 2-9(B) illustrates the substrate for the 
four-point bend coated fatigue specimens. Both coating materials will 
be employed for testing. Table 2- VI lists the number of specimens 
fabricated to date. 
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A) TENSILE, RELAXATION, AND STRESS-RUPTURE SPECIMEN FABRICATED FROM HOT 
ISOSTATICALLY PRESSED POWDER 



B) TENSILE, RELAXATION, AND STRESS-RUPTURE SPECIMEN FABRICATED FROM PLASMA 
SPRAYED SHEETS 



U.w in# \ * L. 

GAGE LENGTH* QIA. 


Figure 2-8 Specimen Designs for Bulk PWA 286 Coating Material Mechanical 
Property Tests 
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A) TENSILE AND CREEP DIFFERENCE METHOD TEST SPECIMEN 



0.25 iron OR 0 J 3 mm 



0.25 cm 
(0.10 in) 


B) FOUR-POINT BEND FATIGUE TEST SPECIMEN 


7.62 cm 
(3.00 in)“ 


1 .27 cm 
(0.50 in) 


l — — 

0.25 


(0.10 


4 

1 

J . 



0.43 mm 
(0.017 in) 


R 


Figure 2-9 Specimen Designs of the Substrates for Coating Material 
Mechanical Property Tests 
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Table 2-VI 


Number of Coating Material Mechanical 
Property Test Specimens Fabricated to Date 


Type of Specimen 

Material 

Number 

Comments 

Cylindrical Tensile, 
Relaxation and Stress- 
Rupture Specimen 

Bulk PWA 286 

29 

Fabricated from 
HIPed powder 

Flat Tensile, Relaxation 
and Stress-Rupture 
Specimen 

Bulk PWA 286 

16 

Fabricated from 
plasma sprayed 
sheets 

Tensile and Creep 
Difference Method 
Specimen 

Bulk PWA 273 

21 

0.25 mm (0.010 
in) substrate 
thickness 


Bulk PWA 273 

17 

0.13 mm (0.005 
in) substrate 
thickness 

Four-Point Fatigue 
Specimen 

Bulk PWA 273 

15 

— 


2.3 Physical, Thermal and Monotonic Mechanical Properties 

2.3.1 PWA 1480 Physical and Thermal Properties 

Physical and thermal property tests of PWA 1480 have been run at Southern 
Research Institute. Data on thermal conductivity, thermal expansion, specific 
heat and density have been obtained and the results are being reviewed. 

2.3.2 Tensile Properties 

2. 3. 2.1 PWA 1480 Single Crystal 

A total of forty monotonic tensile tests were conducted on PWA 1480 single 
crystal specimens with orientations of <100>, <1 1 0>, <1 1 1> , and <123>. All 
tests were run at the ASTM (American Society for Testing Materials) standard 
strain rate of 0.005 min"'. Additional tests are planned at other strain 
rates. Tests included uncoated and aluminide diffusion and NiCoCrAlY overlay 
coated <100> and <1 1 1> oriented specimens. A summary of all test conditions 
and observed material properties is presented in Table 2-VI I . 
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Table 2-VII 

Summary of PWA 1480 Tensile Testing 


Temp. 

Spec. 


Coat 

ExlO” 

. 2 % Yield 

Ult. 

Elong 

RA 

°C( *F) 

ID 

Orient 

Type 

MPa(KSI) 

MPa(KSI) 

MPa(KSI) 

S 

% 

427(800) 

JA-16 

100 

... 

113.8(16.5) 

989.4(143.5) 

1118.4(162.2) 

5.7 

3.2 


KA -2 

110 

— 

221.3(32.1) 

921.9(133.7) 

957.0(138.8) 

14.3 

23.2 


LA-36 

111 

— 

239.3(34.7) 

897.0(130.1) 

1393.5(202.1) 

11.7 

11.9 


MA-1 

123 

— 

198.6(28.8) 

837.7(121.5) 

1218.3(176.7) 

19.1 

18.2 

649(1200) 

JA-33 

100 



POROSITY FAILURE 





KA-3 

110 

— 

176.5(25.6) 

929.4(134.8) 

1081.1(156.8) 

4.7 

5.6 


LA-51 

111 

— 

253.7(36.8) 

849.5(123.2) 

1245.2(180.6) 

23.7 

25.6 


MA-3 

123 

— 

193.7(28.1) 

824.0(119.5) 

1082.5(157.0) 

22.7 

25.1 

760(1400) 

JA-34 

100 

273* 1 * 

286 1,J 

101.4(14.7) 

1177.0(170.7) 

1324.5(192.1) 

14.1 

13.0 


JA-22 


103.4(15.0) 

1159.7(168.2) 

1293.5(187.6) 

4.8 

7.3 


JA-U 


94.5(13.7) 

1163.2(168.7) 

1290.1(187.1) 

8.0 

15.1 


KA-4 

110 


174.4(25.3) 

948.1(137.5) 

1108.7(160.8) 

10.5 

23.9 


LA-52 

111 

273. 1 '* 

200.0(29.0) 

879.8(127.6) 

1093.5(158.6) 

22.1 

24.2 


LA-25 


220.6(32.0) 

920.5(133.5) 

1030.1(149.4) 

16.8 

20.9 


LA-13 


286 i,J 

171.7(24.9) 

908.1(131.7) 

1106.6(160.5) 

21.4 

29.7 


MA-4 

123 

— 

180.0(26.1) 

891.5(129.3) 

985.3(142.9) 

17.8 

18.5 

871(1600) 

JA-36 

100 

— 1 9 

102.0(14.8) 

715.0(103.7) 

1021.1(148.1) 

13.7 

23.6 


JA-24 


273} *‘ 

92.4(13.4) 

756.4(109.7) 

991.5(143.8) 

18.7 

25.4 


JA-12 


286 i,J 

91.7(13.3) 

755.7(109.6) 

961.9(139.5) 

21.7 

26.2 


KA -6 

no 


149.6(21.7) 

786.0(114.0) 

910.1(132.0) 

13.1 

26.0 


LA-53 

111 

273}*? 

190.3(27.6) 

696.4(101.0) 

819.8(118.9) 

19.0 

22.1 


LA-26 


201.3(29.2) 

682.6(99.0) 

812.2(117.8) 

20.3 

24.4 


LA-14 


286 1>J 

181.3(26.3) 

671.6(97.4) 

812.2(117.8) 

22.1 

22.1 


MA-5 

123 

— 

179.3(26.0) 

626.1(90.8) 

764.7(110.9) 

18.0 

21.1 

982(1800) 

JA-37 

100 

273J-2 

286 1,J 

88.3(12.8) 

452.3(65.6) 

695.0(100.8) 

23.0 

32.5 


JA-25 


92.4(13.4) 

437.1(63.4) 

659.9(95.7) 

24.0 

34.3 


JA-13 


102.0(14.8) 

428.9(62.2) 

642.6(93.2) 

22.9 

35.4 


KA -8 

no 

... 

133.1(19.3) 

519.9(75.4) 

628.8(91.2) 

16.7 

36.1 


LA-54 

111 

273*’? 

189.6(27.5) 

427.5(62.0) 

557.8(80.9) 

22.2 

28.0 


LA-28 


175.1(25.4) 

448.2(65.0) 

575.7(83.5) 

18.3 

26.3 


LA- 16 


286 i,J 

120.0(17.4) 

455.1(66.0) 

557.1(80.8) 

23.7 

28.2 


MA -6 

123 

— - 

164.8(23.9) 

431.6(62.6) 

539.9(78.3) 

25.9 

25.0 

1093(2000) 

JA-38 

100 

273}^ 

72.4(10.5) 

275.1(39.9) 

371.6(53.9) 

30.0 

53.8 


JA-26 


52.4(7.6) 

272.4(39.5) 

368.2(53.4) 

31.3 

56.8 


JA-14 


286 i,J 

68.9(10.0) 

269.6(39.1) 

353.7(51.3) 

43.0 

59.2 


KA-9 

no 

— 

91.7(13.3) 

315.8(45.8) 

385.4(55.9) 

18.7 

28.1 


LA-55 

111 

273j*2 

286 i,J 

132.4(19.2) 

259.9(37.7) 

328.9(47.7) 

41.7 

35.5 


LA-29 


97.9(14.2) 

253.0(36.7) 

315.1(45.7) 

28.0 

42.1 


LA-18 


85.5(12.4) 

262.7(38.1) 

321.3(46.6) 

29.2 

46.6 


HA-9 

123 

— 

125.5(18.2) 

273.0(39.6) 

319.2(46.3) 

24.9 

35.3 


1 X-sectional area used to calculate stress excludes coating area 

2 Aluminide diffusion 

3 NiCoCrAlY overlay 
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Figure 2-10 Fracture Surfaces of <001 > PWA 1480 Tensile Specimens. Note the 
pronounced faceting at 760°C (1400°F) is reduced with increased 
temperature. 
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Test specimen JA-33, with orientation <100>, exhibited extremely low 
elongation. Upon metallurgical evaluation, the specimen was found to contain 
large porosity sites along the fracture surface. Such porosity is considered 
noncharacteristic of PWA 1480 material. For actual turbine hardware, the 
standard part inspection practices would detect this type of defect prior to 
installation into an engine. 

2.3.3 Creep Properties 

2. 3. 3.1 PWA 1480 Single Crystal 

To date, eleven of forty planned mo no tonic creep tests have been completed. A 
summary of test conditions and observed material properties is presented in 
Table 2-VIII. 


Table 2-VIII 

Summary of PWA 1480 Creep Testing 


Secondary 


Temp. 

°C ( °F) 

Spec . 
ID 

Orient 

Coat 

Type 

Stress 
MPa (KSI ) 

% Of 

. 2 % Yield 

Life 
(Hrs . ) 

Creep Rate 
(Min' 1 ) 

£1 ong 
% 

RA 

% 

871(1600) 

JA-40 

100 

— 

413.7(60) 

57.8 

462.9 

8.84 E-07 

12.0 

20.1 


JA-41 


— 

517.1(75) 

72.3 

79.0 

8.28 E-06 

15.3 

28.9 


LA-57 

111 

— 

482.6(70) 

69.3 

67.1 

2.90 E-05 

14.9 

22.8 

982(1800) 

JA-42 

100 

--- 

220.6(32) 

48.8 

STOPPED 0 5.4 HRS 

. FOR TEM' 

-1 


JA-27 


273 2,3 

231.7(33.6) 

53.0 

89.1 

3.94 E-06 

25.2 

44.2 


JA-15 


286 2>4 

237.9(34.5) 

55.5 

105.5 

3.65 E-06 

20.0 

42.3 


JA-45 


— 

248.2(36) 

54.9 

80.5 

4.74 E-06 

20.7 

41.4 


JA-28 


273 2,3 

260.6(37.8) 

59.6 

53.3 

8.74 E-06 

24.3 

40.6 


JA-17 


CM 

to 

CO 

CM 

268.2(38.9) 

62.5 

51.7 

8.48 E-06 

20.7 

36.3 


MA-13 

123 

— 

248.2(36) 

57.5 

76.6 

3.12 E-06 

23.5 

33.1 

1093(2000) 

JA-30 

100 

273 2 * 3 

122.7(17.8) 

45.1 

76.4 

2.88 E-06 

20.6 

58.1 


1 Transmission Electron Microscopy 

2 X-sectional area used to calculate stresses excludes coating area 

3 Aluminide diffusion 

4 NiCoCrAlY overlay 
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982°C/220 MPa (1800°F/32 ksi ) to examine the dislocation patterns developed at 

the end of primary creep. At this stage, the dislocation density is low and 



inhomogeous in distribution. Figure 2-11 (A). Dislocation pairs, one is shown 

■ 1 « 1 1 / r\ \ i t ■« .■ m . . I>ie ,1 

the 


ai rs 
f the 


7220 

iqe) 




The <123> oriented creep specimen tested at 'J82°C/248 MPa (1800°F/36 ksi) 
exhibited virtually no primary creep as compared to the <T00> oriented 
specimen (Figure 2-12). Examination of the fracture surface also indicated 
more of a cleavage fracture (i.e., stage I cracking) than is usually observed 
at 982°C ( 1800°F ) . 

A comparison of the percent elongation observed in uncoated <100> oriented 
specimens during rnonotonic tensile and creep testing is presented in Figure 
2-13. As shown, the two elongations are virtually equivalent. 



Figure 2-12 Uncoated PWA 1480 Creep 982°C/220 MPa (1800°F/36 ksi) 
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Figure 2-73 Elongation Observed in Uncoated <700> PWA 7430 Tensile and Creep 
Tests 


2. 3. 3. 2 PWA 286 NiCoCrAlY Overlay Coating 

A summary of the test conditions and observed material properties is presented 
in Table 2-IX. No previous creep experience was available with this material. 
The creep test conditions were set based on the limited stress relaxation 
tests conducted for the constitutive modeling effort. As a result, most tests 
required uploading or were discontinued before rupture. Due to the poor 
success in running these tests and the limited number of bulk hot 
isostatically pressed PWA 286 specimens remaining for constitutive model 
development, the balance of the planned creep tests have been postponed until 
the results of coated PWA 1480 fatigue tests can be used to define new test 
conditions for PWA 286. 
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Table 2-IX 


Summary of Bulk Hot Isostatically Pressed 
PWA 286 Creep Testing 


Temp. 

Spec. 

Stress 

Life 

•C(*F) 

ID 

MPa(KSI) 

(Hrs) 

649(1200) 

9-T 

68.9(10) 

1700 


9-B 

103.4(15) 

1130 

760(1400) 

17-T 

34.5(5) 

446 


12-T 

20.7/55.2(3/8) 

92.1 

871(1600) 

18-B 

6.9/13.8(1/2) 

280.5 


17-B 

20.7(3) 

26.8 

982(1800) 15-B 

3.45( .5) 

— 


Secondary 
Creep Rate 
(Min' 1 ) 

Elong 

% 

RA 

% 

2.99 E-07 

N/A 

N/A 

1.25 E-06 

N/A 

N/A 

6.53 E-06 

93.2 

87.1 

2.64 E-06/4.60 E-05 

166.1 

23.7 


9.38 E-07 /2 .31 E-05 

77.1 

84.9 

1.40 E-04 

206.0 

139.1 

86.7 

66.1 


N/A = Not available 


Comments 
Discontinued 
At 1700 hrs. 
Discontinued 
At 1130 hrs. 

Uploaded from 
20.7 MPa/3 ksi 
to 55.2 MPa/8 
ksi at 48 hrs. 
Uploaded from 
6.9 MPa/1 ksi 
to 13.8 MPa/2 
ksi at 160 hrs 

Failed on 
loading 
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SECTION 3.0 


TASK II - SELECTION OF CANDIDATE LIFE PREDICTION AND CONSTITUTIVE MODELS 


Basic to any life prediction model is a good constitutive model. This is 
particularly true when predicting lives within the hot section of a gas 
turbine engine. As a gas turbine part is cycled through a wide range of 
stresses, strains, and temperatures, deformation and damage accumulate by a 
variety of mechanisms, all of which could play an important role in the part's 
ultimate failure. It is the goal of constitutive modeling to predict this 
stress-strain history and cycle shape, and possibly, the damage history, so 
that the conditions at fatigue initiation are accurately known. For coated 
airfoils, fatigue often originates in the coating and, under certain 
conditions, the coating reduces the fatigue life of the part. Due to this 
intimate interaction between the coating and substrate in fatigue, both 
materials must be understood in terms of constitutive and life models. A 
literature search of numerous such models is reported in this section together 
with the selection of the most promising candidates. 

3.1 Coating Constitutive Models 

Unlike the single crystal substrate materials selected for this program, the 
two coatings, are polycrystalline, and, to a first approximation, isotropic. 
Some anisotropy (Figure 2-4) does exist due to the application process and/or 
interdiffusion between the coating and the substrate. However, for this 
program anisotropy will be ignored for the following reasons: 

1) Coating properties are very difficult to measure and many of the 
planned test specimens are machined from billets which are inherently 
isotropic. 

2) The properties parallel to the surface are believed to be most 
important, while errors in those normal to the surface are 
inconsequential. 

Over this range of complexity five promising models have been selected for 
evaluation in this program. These are: 

1) Classical model, 

2) Walker' s model, 

3) Simplified Walker's model, 

4) Simplified unified approach, and 

5) Stowe 11 equation. 

Each model will be evaluated based on the desirable features which are 
categorized below: 

1) Simplicity 

a) number of model temperature dependent constants 

b) simplicity of test quantity of baseline data required to obtain 
model constants 

c) ease of implementation into a finite element code 
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2) Predictive capability 


a) accurate representation of baseline data 

b) prediction of both isothermal and nonisothermal thermomechanical 
fatigue verification tests 

Priority in ranking the models is depicted graphically in Figure 3-1 together 
with a qualitative preliminary evaluation of each model. Simple models which 
exhibit high predictive capability are considered most desirable and are 
represented by Region I on the figure. Region IV represents the least 
desirable class of models. Emphasis is placed on prediction capability. Thus, 
models which are more complex, such as Walker's unified approach, compare 
favorably to simpler models since the additional complexity should improve 
predictive capability. Quantification of the tradeoffs will become possible as 
test data are generated in Tasks III and V. 


S3 

S3 


REGION III 

REGION 1 

S3 

S3 


m 

REGION IV 

REGION II 

PREDICTIVE 

CAPABILITY ► 


CD CLASSICAL 
S3 WALKER 

S3 SIMPLIFIED WALKER 
S3 SIMPLIFIED UNIFIED APPROACH 

S3 STOWELL EQUATION 


Figure 3-1 Comparison of Constitutive Model Features 


To facilitate quantitative assessment, models are regressed and evaluated 
utilizing a consistent set of isothermal baseline data manipulated by a 
computer automated regression/prediction system. For regressions, the system 
calculates a standard deviation for each model constant which is interpreted 
as a measure of the model's ability to reproduce the baseline data behavioral 
characteristics. 
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3.1.1 Classical Model 


The classical approach (Reference 1) was one of the first attempts at 
developing a nonlinear model which recognizes the observed dissimilarity 
between monotonic tensile and creep inelastic material responses. Time 
independent plasticity (tensile) and time dependent inelasticity (creep) are 
considered as uncoupled components of the total inelastic strain. 

G inelastic = e plastic + E creep ^ 3 


Thus, the total strain function, neglecting thermal strain, is written: 


e total = Elastic + e plastic + £ creep 


(3-2) 


or 


total = o/E + f(o) + 9{ ° st) 0-3) 


where: a = stress, 

t = time, and 
E = elastic modulus. 

Both plastic and creep strain functions [f(a), g(a,t)] are chosen to provide 
adequate duplication of the material behavior. From Task I tests of PWA 286, 
it has been determined that both functions can be described by simple power 
law relationships: 


/ \ c 

I CT 1 

e plastic = \K^ I 

A 4 

E creep = t 


where: Aj, A 3 , and A 4 are constants. 

These forms will be used in the first curve fit of constitutive data and will 
be modified as necessary as more data become available. 
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3.1.2 Walker Model 

The Walker model (Reference 2) is among a new generation of constitutive 
models based on a unified viscoplastic approach which considers all nonlinear 
behavior as time dependent inelasticity. No distinction is made between 
plastic and creep inelastic action as in the classical model. Walker, from his 
earlier work on Haste! loy X, has chosen to express inelastic behavior by a 
power law relationship which can be written one-dimensionally as: 



finelastic 


\-e~l 


(3-6) 


where n is a constant and £2, back stress, and K, drag stress, are strain 
history dependent internal state variables which describe kinematic and 
isotropic cyclic hardening, respectively. 

The back stress term, Q, is a deviatoric quantity which physically corresponds 
to the asymptotic stress state under relaxation conditions. Qualitatively, the 
evolutionary expression for back stress is a sum of opposing nonlinear 
hardening and thermal and dynamic recovery components which can be 
characterized as: 


a = f (e in , e in , T, t) - f (O, e. n , T, t) (3-7) 

Hardening Recovery 

where: e,j n = instantaneous inelastic strain, 

£i n = instantaneous inelastic strain rate, 

T = temperature, and 

t = time. 

Drag stress is a scalar quantity which represents a resistance to inelastic 
flow, and is considered a function of the effective inelastic strain, R. 

One-dimensionally: K = K 1 ' K 2 exp (“ n 7~) (3-8) 

where: K * instantaneous drag stress, 

= fully hardened/softened drag stress, 

«i - l <2 = initial drag stress, 
ny = constant, and 
n 

R = .2^ |Ae^ n I ; Ae^ p = inelastic strain increment. 
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Thus, the drag stress function is a monotonically increasing relationship 
describing isotropic hardening (K 2 > 0) or softening (K 2 < 0). 

The entire set of equations for the Walker model considered in this research 
effort are provided in Table 3-1. References 2 and 3 provide a detailed 
discussion of Walker's and other unified approaches. 


Table 3-1 

A Summary of Walker's Model 


• / V2/3 (3/2 s..- ,.) (3/2 s.j-OtjA " 


irV 


y 2 / 3 (3/2 s ir n..) (3/2 s^j) 


°ij = ^ij * c 1j TiT 8 * n l®fj * *\ l i* "ij 1 


Vij- °ij l) («i -ig-S °) 

■ -u - "Sr (i 2 - 4 ^ ;) 

V ("3 * •*(*«!}’ “ij") 

=2 ■ "9» * "10 (f "if “if ) 


H-l 


-n 7 R -n„L 
K = K : - K 2 e ' - K 3 e 0 


» . 1 / 2/3 

2 ■ K 2/3 "l3'lj - I ^ ‘ij "IJ I 


s ij * °1j * T 40 kk 

"--‘ij s ij* 34 ucjjc*J 


'Ij 


s lj x u kk 


• • • * 
c kk + 2 " c 1j - °1j * s ij (3X + 2 “ ) J “ 8 d 


/ 2u 


Material Constants: X, u, n, n, m, ni, n 2 , n 3 , n 4 , ns, ng, n 7 , 
H8« n 9» ni0» n ll» K l» k 2» an d k 3 depend on 
temperature. 


( 1 ) 

( 2 ) 

(3) 

( 4 ) 

(5) 

( 6 ) 

(7) 

( 8 ) 

(9) 

( 10 ) 

( 11 ) 

(12) 
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3.1.3 Simplified Walker Model 


Preliminary Task II results suggest that the Walker model can be simplified 
because some terms describing the evolutionary response of back stress are 
unnecessary for PWA 286 coating material. The simplified Walker model is an 
extension of this simplification process to the extreme case where all back 
stress terms are eliminated. This seems justified, at least qualitatively, 
based on the rapid relaxation/low strength properties of coating materials at 
moderate temperature levels, T _> 871 U C (1600°F). Thus, any exposure at high 
temperature can erase the back stress component (Figure 3-2). The inelastic 
strain rate formulation, equation (3-6), is then reduced to: 



(3-9) 


From a simplicity standpoint this model is very attractive. 



TEMPERATURE *C (*F) 


^ CALCULATED BY WALKER MODEL SIMULATION OF A MONOTONIC 
TENSILE TEST: e= 2%; e= 0.005 min"' 


Figure 3-2 Saturated Back Stress Levels for Unexposed PWA 286 Overlay 
Coating 
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3.1.4 Simplified Unified Approach 

Recently, Moreno (Reference 4) has had success in predicting the constitutive 
response of Hastelloy X using a hybrid model derived from both classical and 
unified approaches. Basically, the model assumes that the thermo-elastic 
response is known from previous analysis from which incremental values of 
strain, temperature, and time are used to calculate the actual (nonlinear) 
stress history. 

Strain increments are considered as either time independent plastic or 
thermo-elastic creep. 


Ae = Ae . . 

plastic 


(3-10A) 


or 


Ae A Elastic + " e creep 


(3-10B) 


Expressed in terms of stress increments: 


Ao = Act 


plastic 


(3-11A) 


or 


Act = ACT elastl . c + creep (3-11B] 

The classical yield surface criteria (Figure 3-3) are used to determine the 
onset of plastic action. Isothermal yield points are calculated from a 
tri-1 inear representation of tensile stress/strain information (Figure 3-4), 
and are assumed to be equal in tension and compression (i.e., the material is 
perfectly isotropic). No cyclic hardening is considered. Justification for 
this definition of yielding was based on observations of Hastelloy X cyclic 
response: 1) at elevated temperatures, little cyclic hardening/softening 
occurs, and 2) during thermomechanical cycling, exposure to the higher 
temperatures significantly reduces the amount of cyclic hardening developed at 
the lower temperatures. 

The stress increment associated with time independent plastic action (equation 
3-11A) is then calculated as: 

Vl " a i = Ao plastic = a y. +1 " a y. (3-12A) 


for 


a . = a 

i 


T ,n> T , 
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or 


E + E 

PiH Pi 

°i +1 " CT i ~ Aa plastic = T “ Ae 


for 


a . = cr 


T i+t < T i 


(3-12B) 


and the stress increment associated with time dependent creep behavior 
(equation 3-11B) is calculated as: 


i + 1 


= Act 


elastic + creep 


= E ( Ae 


A e ) 

creep 


(3-13) 


A e 

creep 


Aa n At 


(3-14) 


where: A(J elastic + creep = total thermo-elastic creep stress increment, 

Aapi d stic = total time independent plastic stress increment, 

0 y = yield stress, 

Ae = total strain increment, 

Ae C reep = total creep strain increment, 

E = elastic modulus, 

Ep = strain hardening slope of monotonic a/e curve, 

T = temperature, 
t = time increment, 

subscripts: i = beginning of increment, and 

i+1 = end of increment. 

This model has been successful with Hastelloy X, which exhibits some 
behavioral similarities with overlay coating PWA 286. However, the higher 
level isothermal cyclic hardening characteristic of PWA 286 may be more 
difficult to simulate with this simplified model. 
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STRESS 



Figure 3-3 Classical Temperature Dependent Yield Surface for Simplified 
Unified Approach 



STRAIN 

Figure 3-4 Construction of Tri-Linear Stress Strain Curve 
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3.1.5 Stowell Model 


The Stowell model (References 5 through 7) is another form of a unified 
viscoplastic approach. It considers inelastic action as a diffusion mechanism 
based upon an apparent activation energy level. 


£ inelastic 


= 2S T exp 



(3-15) 


where: 


aH = apparent activation energy, 

a = applied stress, 

T = absolute temperature (°K or °R), 

R = gas constant, and 

S, a 0 = constants. 


Since creep deformation is generally thought to be controlled by diffusion 
processes, this model should perform well at elevated temperatures and/or slow 
strain rates. 


3.2 Single Crystal Airfoil Material Constitutive Model 

In Task II, a literature survey (Appendix B) was conducted to identify 
constitutive models and modeling approaches that are applicable to single 
crystal and transversely isotropic materials. Two distinct modeling approaches 
can be identified. The first attempts to model the anisotropic behavior by 
using the micromechanical processes which occur during deformation as a basis 
for the mathematical formulation. The second approach ignores the 
micromechanical processes and attempts to describe the bulk behavior through 
macroscopic quantities which are made to be functions of material orientation. 
While some of the earlier macroscopic models can be considered "fully 
developed", no complete micromechanical models can be considered fully 
developed. 

Three candidate models have been identified for further evaluation. Two of 
these models use the macroscopic continuum approaches. They are the classical 
Hill model which is currently coded in general purpose finite element codes 
(such as MARC and ABAQUS), and a unified viscoplastic formulation of a model 
presented by Lee and Zaverl , et al . The third model uses a micromechanical 
approach which is currently under development. Each of these models are 
discussed in the following sections. 

In interpreting the test results extensive use is made of the stereographic 
projection to display the angular relationships between crystallographic 
directions. All the planes in a crystal can be represented by a set of plane 
normals radiating from a single point of origin within the unit cell of the 
crystal. If a reference sphere is described about this point of origin, the 
plane normals can be described as the points where the normals intersect the 
surface of the sphere. This process is illustrated in Figure 3-5 for just the 
jlOO} planes of the cubic unit cell. 
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Figure 3-5 . The jlOO} Poles for the Cubic System 


The stereographic projection is illustrated in Figure 3-6 which shows a small 
cube of face centered material located at the center of a sphere. The axes CC 
is aligned in the Z or [001] direction, while the axes BB and AA are aligned 
in the y or [010] direction and in the x or [100] direction, respectively. 
Figure 3-7 shows the orientation ON which is defined by the Euler angles e and 
^ with respect to the cubic crystallographic axes x, y, z. If this direction 
ON is projected onto the sphere in Figure 3-6, it intersects the spherical 
surface at point T. The intersection of the line CT with the diametral circle 
ABAB occurs at point E, the circle ABAB being the stereographic projection 
with the [001] axes represented as the line OC. The direction OL where e = 45° 
and ip= 0° in Figure 3-6 intersects the sphere at point P and the line CP, 
which goes through point L on the corner of the cube, intersects the 
stereographic projection at point R. This i$ o the [Oil] direction. Similarly, 
the direction OM where e = 45° and >P= 35.26° intersects the sphere at point Q 
and the line CQ intersects the stereographic project at point S. This is the 
[111] direction. The stereographic triangle ORS thus represents all angular 
orientations with respect to the crystallographic cubic axes. From the 
symmetry of the cube any other direction outside of the foregoing limits is 
equivalent to a corresponding direction inside the limits of the stereographic 
triangle ORS. 
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Figure 3-6 A Single Bar Oriented Along the Direction ON With Respect to the 
Cubic Crystal Axes AA, BB, and CC can be Located in the 
Stereographic Triangle ORS at Point E. The circular area ABAB is 
the stereographic projection. 
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mu 

it 45 *.* = 35 . 26*1 



Figure 3-7 A Single Crystal Bar Oriented Along ON at Angles e and ip With 

Respect to the Crystal Axes X, Y, Z is Located at Point E in the 
Stereographic Triangle ORS. 
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3.2.1 Classical Hill Model 


In 1948, Hill (Reference 8) proposed a generalization of the von Mises yield 
function for materials having three mutually orthogonal planes of symmetry. 
This yield function and an associated flow rule are the most widely used 
formulation for anisotropic materials and has been coded into general purpose 
nonlinear finite element codes such as MARC and ABAQUS. In this sense, it 
represents "standard practice" for analyzing structures with anisotropic 
material. For this reason, it was selected for further evaluation in this 
contract. 

The Hill yield function is most commonly written in the form: 


- I3\ 1/2 I 1 \ 1/2 


F( CTy - CT Z ) 2 + G(o z - a x ) 2 + H(o x - 0y) 2 


(3-16) 


+ 2L T 2 + 2M r 2 + 2N t 2 
yz zx xy 


1/2 


where a is the effective stress and a x , ay, a z , etc. are the stresses 
referred to the material coordinate system. F, G, H, L, M, and N are constants 
which in effect permit yielding to occur at different stresses in the 
different material directions. For single crystals, the material can be 
assumed to behave identically along the three crystal axes, [100], [010], 
[001], so that symmetry arguments can be used to show that: 


F = G = H and L = M = N. 
The effective stress can then be written as: 

'3 \ 1/2 II \ 1/2 




F|( "y - °z >2 * <°z - °x )Z * ,0 x 


“ 0 , 


■> 2 } 


h 2l It 2 + r 2 + r 2 } 
l yz zx xy/ 


1/2 


(3-17) 


or 


2 + 2 + 2 
°x CT y a z 


aa - a a - a a + A 

y z z x xy 


2 . 2 . 2 
r t + t 
yz zx xy - 


1/2 


(3-18) 


where A = L/F. 
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For the case of a simple uniaxial specimen oriented arbitrarily In the 
crystal, the stresses referred to the crystal axes are: 


2 2 

o = cos ip cos e S 

A 

2 2 

0 = cos 4 1 cos e S 
a z = cos 2 ^ S 


2 

7 xy = cos 4> cos e sin e S 

T yz = cos2 ^ sin ip sin © S (3-19) 

r = cos tp sin >P cos © S 

bA 


where e and */> are the Eulerian angles shown in Figure 3-7 and S is the 
uniaxial stress in the specimen. 

Substituting, the effective stress equation becomes: 

0 = S jf(^, ©) - g {ip, e) + l\g(ip, ©) } 1/2 (3-20) 

where 


f(ip, ©) = cos 4 'Z' (cos 4 © + sin 4 e) + sin 4 '/' 

(3-21) 

and g (ip, ©) = cos 4 cos 2 © sin 2 © + cos 2 '/' sin 2 '/' 

Equation (3-20) defines the orientation dependence of the yield stress 
permitted by the Hill model. For a uniaxial specimen oriented in the [001] 
direction (</'=© = 0), equation (3-20) gives the effective stress at yield: 


a = S [001] 


(3-22) 


For a specimen oriented in the cube diagonal direction [111] (© = 45°, if/ = 
35.26°), the effective stress at yield is: 


a = 



(3-23) 


*[ 100 ] 

5 [ 111 ] 


A = 3 


2 


(3-24) 



Substituting equation (3-24) into (3-20) and manipulating gives an expression 
for the yield strength in any direction S as a function of the yield stress in 
the [001] and [111] directions. 


S , S [001] | fit ') - git o) * 3 9lt 9) | (3-25) 

Figure 3-8 shows the permitted variation in yield stress for the hypothetical 
case where the cube diagonal yield stress is twice the cube axes yield stress. 
This figure shows that the permitted variation in yield strength changes 
monotonically from the [001] to the [111] direction which is not in agreement 
with reported single crystal data. This represents a severe restriction in 
this classical model. 

Further limitations of the Hill formulation have also been identified. The 
yield function is insensitive to the sign of the stress state so that it 
cannot properly account for observed differences in tension and compression 
data (Reference 9). In addition, the state of anisotropy does not change with 
deformation. In spite of these obvious limitations, the classical Hill model 
will continue to be evaluated against the body of data being obtained because 
of the model's widespread use. 



Figure 3-8 The Hill Yield Function Permits Only a Monotonically Changing 
Yield Strength in the Stereographic Triangle 
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3.2.2 Macroscopic Viscoplastic Formulation 

The second modeling approach selected for evaluation and development is a 
generalization of the rate- independent yield function proposed by Lee and 
Zaverl, et al (Reference 10). For application to structural analysis of 
turbine airfoils, a rate- independent formulation is not as desirable as a 
unified viscoplastic formulation. Therefore, this second model will imbed the 
anisotropic matrices D-jj and Mjj of Lee and Zaverl in a unified 
viscoplastic framework. 

Isotropic unified theories were reviewed in a pre-HOST program in Reference 
11; these theories are readily modified to include anisotropy. For example, 
the unified theories may be modified so that the expression for the inelastic 
strain rate, Cjj, has the form: 


ijjkl (a ij- n ij } (a kl" n kl } 
r— ■ =ZI~ (3-26) 

v i'jkl {o ij" J2 ij ) (a kT n kl ) 

where J2jj and K are the equilibrium (back or rest) stress and the drag 
stress, respectively. The power law expression in equation (3-26) may be 
replaced by an exponential, hyperbolic sine, or any suitable functional 
relation according to the particular theory under consideration. The matrix 
form M-jj of the fourth rank tensor Mij^i has two independent components 
for single crystal materials and three independent components in the case of 
directionally solidified materials. 

The growth law governing the evolution of the equilibrium stress, 12^,-, may 
be assumed to have the form: J 


VmTj 


c ij = 


ijkl (<y ij~ n i.j ) (g kl~ fi kl ) ^ 


where 


n ij = a l c ij “ (n ij 


fl,- A “ ® * 


1j 


1J 


h = the initial value of the back stress 
G = a 3 R a 4 (M 1Jkl fi.j. S2 kl ) ^r 




( a i j " n ij^ (a kl " fi kl ^ 
K 



(3-27) 

(3-28) 

(3-29) 


The drag stress, K, may be assumed to evolve according to the equation: 


K = a 5 R - a g KR - a ? K q , 


(3-30) 


with R defined in equation (3-29). 
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Calculations may be made with this anisotropic unified theory to construct a 
theoretical yield surface. For a given point in stress space, a surface 
surrounding this point may be constructed such that in going from the given 
point to each point on the surface, the cumulative inelastic strain increment 
is the same for each point on the surface. If this theoretical construct is 
called the yield surface, then the model allows the yield surface to translate 
in stress space by virtue of the presence of the equilibrium stress (similar 
to kinematic hardening in the isotropic version) and to expand in stress space 
due to the presence of the drag stress (similar to isotropic hardening in the 
isotropic version). In addition, the yield surface can distort its shape and 
rotate in stress space due to the presence of the fourth rank anisotropy 
tensor M-ijki which might be assumed in the form: 



= (M, 


ijkl 


" M jkl } 


(3-31) 


where M ijkl denotes the limiting saturated,, value of the tensor as 

r— oo under cumulative deformation, with M ijkl = Mijkl being the initial 
value of Mijki. A thermal recovery term may be added to equation (3-31) at 
high temperatures to allow recovery of the material to its original 
anisotropic state, viz. j^i . 

The unified formulations are completed by the specification of the relation: 


ij 


= D ijkl 


"kl 


c kl> 


(3-32) 


where Djjj,-) is the fourth rank anisotropic elasticity tensor. 


3.2.3 Micromechanical Viscoplastic Formulation 

The third constitutive model selected for further development is a 
micromechanical model based on slip system stresses and strains and imbedded 
in a unified viscoplastic formulation. Development of a slip system based 
constitutive model began at Pratt & Whitney in 1982 and continues to be part 
of a large company effort to develop constitutive models. The slip system 
model development has continued since 1982 through the combined efforts of 
Pratt & Whitney's personnel and consultants under the company program and in a 
complementary NASA grant NAG-512. 

The viscoplastic consitutive formulation based on crystallographic slip theory 
can be written in the form (Reference 12): 
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For Octahedral Slip 

m 1 =( i - k)/vT, m 2 = ( - i + jl/vT, m 3 =(-j + k)/vT, m 4 =(j - k)/ V2 , 
m g =(-i - ji/VT, m 6 =(i + k)/vT, m 7 = ( - i - k)/V2~, m 8 =(i - j)/vT\ 
m 9 =(j + k )/VT, m 10 =(-j - k)/y/2~, m n =(i + j)/VT, m 12 =(-i + k)/\ZT , (3.33) 


n i = n 2 = n 3 = (i + j + k)/y/Ti n 4 = n 5 = n g = (-i + j + k)/ VJ”, 

n 7 = n g = n 9 = (-i - j + k )/VT, n 1Q = n n = n 12 = (i - j + k)/V r 3 _ , 

where X> J > and k are the unit vectors in the Cartesian coordinate system 
having axes parallel to the three coordinate axes. 

= m ‘ c ‘ n, , tL, = ^ i = n / c ‘ n, , 


mn r 


jr ’ mm _r 


~r 7 nn „r 


r . r . . r 

7T = Z„ ’ a * Z„ , ir_, = m * a * Z„ , ir _ = n„ ’ a ’ Z„ , 

zz „r _ jr mz „r „ ~r nz „r „r 


(3-34) 


fr “ !?r ' Ur > 


(3-35) 

r , r . r . 

*r = "mn a mm "mm a nn *nn 

“zz ’zz + 2 “mz ’mz * 2 “nz ’Hz ■ 

(3-36) 

Y r = K~ p ( r, . u r ) | ir r - u r | 

P-1 

3 

(3-37) 

00 • 

“ “ P 1 V ‘ P 2 | Y r | - P 3 

1 1 m_1 
l u r 1 u r * 

(3-38) 

K r = K«*P 4 tanh ('“) * £ 

3K r 

~w * 

(3-39) 

*r = r=l h rk 1 Y kl ’ 


(3-40) 
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(3-41) 


h 


rk = 


0 


q + (i - q) « rk 


e - h 





dS , 


c 


12 1 • 

k!l 7 Y k <Jk ”k * ”k Jk> > 


(3-42) 


o • • 

o = D : (e - c) . (3-43) 

m mm 

For Cube Slip 

uii = (i + m ? = (-1 + j)/vF, m~ = (i + k )/\/2“ , 

m 4 > (-i + k)/vT, ra q = (j + kJ/V^FT m fi = {-j + k)/vT, (3-44) 

MM M^ MM mD MM 

!!l ■ Iz - !;• "3 = ?4 = b ?5 ■ ?6 - i • 


The remaining equations are the same as for the octahedral slip except that 
the summations involving the indices r and 1 range from 1 to 6 instead of from 
1 to 12. 

A preliminary evaluation of this model, has been made using the 982°C (1800°F) 
data obtained to date. Figure 3-9 shows the effect of orientation on Young's 
modulus for PWA 1480 at 982°C (1800°F). The modulus increases from a minimum 
value of 78.6 x 10 3 MPa (11.4 x 10® psi) at the [001] corner of the 
stereographic triangle and reaches a maximum value of 219. 3 x 10 3 MPa (31.8 
x 10® psi) at the [111] corner of the triangle. The elasticity matrix for a 
material which exhibits cubic symmetry can be written in the form: 


11 

°12 

°12 

0 

0 

0 


12 

D 11 

°12 

0 

0 

0 


12 

°12 

D 11 

0 

0 

0 

(3-45) 


0 

0 

°44 

0 

0 



0 

0 

0 

°44 

0 



0 

0 

0 

0 

°44 
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where D n . 206.8 x 10 3 MPa (30 x 10 6 psi), 

Dj2 = 151.7 x 10 3 MPa (22 x 10 6 psi), and 
D44 = 85.5 x 10 3 MPa (12.4 x 10 6 psi) 
for PWA 1480 at 982°C (1800°F). 



Figure 3-9 Stereographic Contour Plot of Young's Modulus for PWA 1480 at 
982 C (1800°F) 


To determine the material constants for octahedral slip, the COPES/CONMIN 
(Control Program for Engineering Synthesis/Constrained Minimization) code 
(Reference 13) was used. The objective function to be minimized in the 
COPES/CONMIN code was chosen to be the main square of the difference between 
the experimental and theoretically predicted cyclic hysteresis loops carried 
out at strain rates of 10 -3 , 10“ 4 and 10 -5 per second in the [111] 
direction. According to the COPES/CONMIN code, the minimum mean square 
deviation between the experimental and predicted hysteresis loops occurs for 
the following values of the material constants: 
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KP = 28317 for r=l to 12, 

P1 = 2.49 x 10 7 , 

P2 = 2871, 

P3 = IQ" 8 , 
p = 6.34, 
m = 2.12, 

P4 = 940, 

P5 = 10400, 

0 = 0 , 

q = 1 (no latent hardening assumed for the present), 
h = arbitrary (since 0=0), 

D n . 30 x 10 6 , 

D 12 = 22 x 10<S, 

D44 = 12.4 x 10 6 , and 

a mm = a nn = a zz = a mz = a nz = 0* 

Figures 3-10 and 3-11 show the experimental and predicted hysteresis loops in 
the [111] direction for strain rates of 10- 3 , 10 -4 and 10 -5 per second. 

For octahedral slip, the variation of the constitutive response with respect 
to direction in the stereographic triangle is now fixed by the preceding 
values of the material constants. Figures 3-12 and 3-13 show the comparison 
between the experimental and theoretical cycle hysteresis loops in the [001], 
[Oil] and [111] directions at a strain rate of 10' 3 per second. The 
variation of peak stress response with direction is reasonable, but the 
theoretical hysteresis loop in the [001] is "fatter" (i.e., displays a more 
inelastic response) than the corresponding experimental loop. At this strain 
range, viz. ^0.3 percent, the experimental [001] loop is almost behaving 
elastically. At higher strain ranges, the response exhibits an inelastic 
behavior and the loop becomes "fatter". 
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e- 10"3 sec"] 
e= 10-]] sec" 
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Strain, percent 


Figure 3-10 Experimental Hysteresis Loops in the [111] Direction for PWA 
1480 at 982 °C (1800°F) 



«= 10“3 sec -1 
f» 10" 4 sec" 1 
e= 10" 5 sec -1 


Figure 3-11 Precited Hysteresis Loops in the [111] Direction for PWA 1480 at 
982°C ( 1800°F ) 


Stress, psl 





Figure 3-14 shows the corresponding predicted loops at a 10"^ per second 
strain rate in the [001], [Oil] and [111] directions using the cube slip 
formula. The same constants were used here for the cube slip response as were 
used for the octahedral response, simply to see the variation of cube slip 
constitutive response with direction. The material constants are clearly 
inappropriate but suffice to show that no inelastic response is predicted in 
the L001 ] direction. Even for large strain ranges, no inelastic response is 
predicted in the [001] direction due to the fact that the resolved shear 
stress r mn is zero in each of the six slip directions when the bar specimen 
is pulled in the [001] direction. Clearly a mixture of octahedral and cube 
slip will improve the correlation between the theoretical model and the 
experimental results shown in Figures 3-12 and 3-13. 



Figure 3-14 Cube Slip Predictions of Hysteresis Loops for PWA 1480 at 982°C 
( 1800°F) 


The octahedral slip formulation was used to predict the response during a 
tensile test at 982°C (1800°F) and 10“3 sec”* strain rate. The results are 
presented in Figures 3-15 through 3-17. 
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Theta, degrees 


Figure 3-17 Stereographic Contours of Inelastic Strain at 50 ksi Stress 


The stereographic contours of the tv/o- tenths offset strain are shown in Figure 
3-15, while Figure 3-16 depicts the stress contours at the large offset strain 
value of two percent. The two percent offset stress contours have an 
appearance which is similar to the variation of the octahedral Schmid stress 
over the [001], [Oil] and [111] stereographic triangle. When the inelastic 
strain contours corresponding to a load of 50,000 psi on the uniaxial bar 
specimen are plotted out in Figure 3-17, it is evident that they follow the 
Schmid stress variation exactly. This is to be expected since an octahedral 
slip theory is being used. However, at low values of the offset strain the 
contours do not follow the Schmid pattern since the state variables u r and 
K r are in a growth phase. In Figure 3-18, the stress contours corresponding 
to a two- tenths offset strain are plotted for the cube slip model where it is 
evident that the stress values decrease from the [001] vertex, where the 
"yield" stress is infinite, towards the [111] vertex. 
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Further work will address the constitutive modeling of the 427°C, 760°C and 
871°C (800°F, 1400°F and 1600°F) data using a mixture of octahedral and cube 
slip. Work is progressing under NASA grant NAG-512 to determine if the 

material constants can be obtained without the use of the COPES/CONMIN code. A 

symbolic computer algebra program called MUMATH, which operates on the IBM 
PC-X and Apple personal computers, is being employed for this work. If the 
material constants can be determined by this approach, the COPES/CONMIN code 
can be used to "fine-tune" these material constants. At present, the constants 
are simply guessed and fed to the COPES/CONMIN optimization code, which 
operates on the IBM PC-XT computer. If the guesses are not close to the 

optimum values, the optimization code will use a large amount of computer time 

in determining the material constants. 



Figure 3-18 


Stereographic Contours of Cube Slip Stress Corresponding to 
Two-Tenths Offset Strain 
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Figure 4-3 Short Time Creep 








Table 4-1 


Summary of Baseline Stress Relaxation Tests for PWA 286 
NiCoCrAlY+Si+Hf Overlay Coating 


Exposed =100 hours at 1093°C (2000°F) 


Test Virgin (V) Temperature 

Number or Exposed (E) *C (°F) 


Low Plasma 

HIP PWA 286 Spray PWA 286 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


vi 

vi 


E 

V 

E 


Vl 


538 (1000) X 2 
760 (1400) X 
760 (1400) X 
871 (1600) X 
871 (1600) X 
982 (1800) X 
982 (1800) X 
1038 (1900) X 
1038 (1900) 

1093 (2000) X 


X 

X 

X 

X 

X 

X 


lllsed for constitutive model evaluations 
2 Each "X" represents a single test specimen 


Table 4-II 

Summary of Baseline Constitutive Tests for PWA 273 
Aluminide Diffusion Coating 

Exposed =100 hours at 1093“C (1900°F) 


Test 

Number 


Virgin (V) Temperature 

or Exposed (E) a C (!f) 


1 

2 

3 

4 

5 

6 

7 

8 
9 


vi 

V 1 

vl 

E 

Vl 


Vl 

E 


538 (1000) 
649 (1200) 
760 (1400) 
760 (1400) 
871 (1600) 
871 (1600) 
982 (1800) 
1093 (2000) 
1093 (2000) 


lllsed for constitutive model evaluations 

The tests will consist of either monotonic or cyclic waveforms. Specimen 
design is shown in Figure 2-9(A). 
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Model constants were determined for the unexposed PWA 286 overlay coating from 
stress relaxation tests run at 538, 760, 871 , 982, and 1093°C (1000, 1400, 
1600,^1800, and 2Q00“F) [see Appendix A). For the model discussions, 538°C 
(1000“F) and 982 J C (1800 F) data are presented as being representative of the 
low and high temperature material response. PWA 273 diffusion coating baseline 
constitutive tests are pending. 

4.1.1 Classical Model 


A 2 A 4 

•ft* (y * (%) 1 

where: et = total strain 

e e = elastic strain 

£p = plastic (time independent inelastic) strain 
e c = creep (time dependent inelastic) strain 
a a stress 
E = elastic modulus 
t = time 

A l» A 2» a 3» and A 4 are constants. 

Baseline data regressions for unexposed PWA 286 utilizing the classical model 
given above are presented in Figure 4-4. As expected, time independent 
plasticity dominates the inelastic response at lower temperatures, while creep 
dominates at higher temperatures. Figure 4-4 also indicates that some history 
dependence is required to accurately represent the observed isothermal cyclic 
hardening. 
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4.1.2 Walker Model 

The one-dimensional form used to regress the baseline data was: 


e t = 

£ e 

+ e in 

• 

1 

a - 0 \ n 

e in 

■ V 

-t~) 

K = 

K i 

- Kg exp ( 

S2 = 

n i 

e in + n 2 


n 2 = n ll e in " a 2 (G 2 


'll 


an 

ae 


©) 



V 


"10 n 2 


<Vl> 


R = 



where: ein = inelastic strain 

e-j n = inelastic strain rate 

n = back stress 

K = drag stress 

R = effective inelastic strain 

R = effective inelastic strain rate 

n, ni, ny, ng, nio, nn, Kj, K 2 , and m 0 are constants. 

The two back stress terms (^i and fig) of th e Walker model (Table 3-1) are 
mathematically similar; therefore, it was felt that one of the terms could be 
eliminated. A study was conducted to determine the regression capability using 
% and £22 together and separately. The results, presented in Figures 4-5 
through 4-7, indicate that no regression capability is lost if either one of 
the back stress terms were dropped. This is also indicated by the standard 
deviations of the regressions which are provided in the lower right-hand 
corner of each figure. Since the term contains more constants than the 
^2 term, that term was eliminated. 
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Figure 4-5 Walker Model Regression Fit of Unexposed PWA 286 Stress 
Relaxation Data Using and Back Stress Terms 
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Figure 4-6 Walker Model Regression Fit of Unexposed PWA 286 Stress 
Relaxation Data Using Back Stress Terms Only 
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Figure 4-7 Walker Model Regression Fit of Unexposed PWA 286 Stress 
Relaxation Data Using Back Stress Terms Only 
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4.1.3 Simplified Walker Model 


e t = 



E 


in 



K = - K 2 exp (-n 7 R) 


R = 



The above term descriptions are the same as for the Walker model (see Section 
4.1.2). 


Figure 4-8 presents the regressions using the Simplified Walker model. 
Clearly, no gross reduction of regression capability is produced by ignoring 
back stress. 





4.1.4 Simplified Unified Approach 


For a < a . , 


e t = f + A a B t 


For a. = a yi . , 


A a 


plastic = a yi+l ' CT yi ; for T i + 1 ^ T i 


where: 


subscripts: 


Act 


plastic 


; Er1 “ t ; for T m < T i 


Aopiastic = time independent plastic stress increment 
CTy = yield stress 

E p = strain hardening slope of monotonic a/e curve 
Aet = total strain increment 
A, B are constants 
i = beginning of increment 
i + 1 = end of increment 


The Simplified Unified Approach model has not yet been regressed because it is 
felt that the model cannot adequately represent the baseline isothermal tests. 
Unlike Hastelloy X, significant isothermal cyclic hardening/softening was 
observed in PWA 286. Regression of this model is pending evaluation of 
nonisothermal tests. If those tests indicate that the net cyclic 
hardening/softening is insignificant, the model will be re-evaluated. 

4.1.5 Stowe! 1 Model 


e 


t = 



e 


in 


e- n = 2sT exp (^y) sinh (—■) 

1 1 ,,, o 
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where: aH = apparent activation energy 

T = absolute temperature (°K or U R) 

R = gas constant (2 cal/mol e. °K or 1545 ft.lb/lb m -mole. ! r) 
S, o 0 are constants 


The baseline regression data fits using the Stowell model are presented in 
Figure 4-9. Note that, at 538 J C (1000°F), the model is unable to duplicate the 
observed yielding during loading. This suggests that yielding at low 
temperature is a nondiffusion process which the Stowell model is not well 
formulated to correlate. 


xl05 


533»C (1000°F) 


w (0.8) 
* ( 0.6 
| (0.4) 

( 0 . 2 ) 
0.0 
(- 0 . 2 ) 
-275.6 (-0.4) 































7 

7 




• 







A 

.. 

7 

/• 

7 


> 









—A 



7 

a y 

V 


. 









A 


•y 

V 

> 

7 


• 

7 








A 

Y 


7 

-V 

7 



7 








A 

T 

V 

7 

7 

/ 



y 








A 

T 

A 

V 

7 



a 

7 








A 

7/ 

A 


7 

7 



7 
















7 




























0.000 


0.002 0.004 0.006 0.008 0.010 0.012 0.014 0.016 

Strain, In/in 

982°C (1800'F) 


Data 

Regression Fit 



Figure 4-9 Stowell Model Regression Fit of Unexposed PWA 286 Stress 
Relaxation Data 
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4.1.6 Summary 


Each model considered was regressed on a consistent basis utilizing a computer 
automated regression system. The standard deviation between the regressed 
model and the data is used to quantify the regression capability of the model. 
A constitutive model summary for unexposed PWA 286 is presented in Table 
4-III. The results indicate that the various models essentially compare as 
presented in Figure 3-1. 


Table 4-III 

Summary of Unexposed PWA 286 Coating Constitutive Models 


Finite Element Code 1 Standard Deviation of Baseline Data Regressions 2 




Similar to 


Temperature °C (°F) 



Model 

Speed Availability 

Substrate 

538 (1000) 

760 (1400) 

871 (1600) 

982 (1800) 

1093 (2000) 

Classical 

X 


3147.0 

1866.0 

748.5 

310.0 

128.0 

Walker 

X 

X 

1999.0 

1816.0 

633.4 

152.8 

100.7 

Simplified 

Walker 

X X 

X 

2041.0 

1717.0 

876.1 

220.4 

118.5 

Simplified 

Unified 

Approach 

X 


• 

“ 

* 

• 

* 

Stowe 11 X 

| Anticipated finite element code qualities 
2 lo; expressed In psl 

6521.0 

2091.0 

1044.0 

376.9 

139.5 


4.1.7 Future Work 

In the remaining period of this task, additional testing, both isothermal and 
nonisothermal , will be conducted to acquire information for: 

1. verification of PWA 286 coating constitutive models, and 

2. regression and verification of PWA 273 coating constitutive models, 
and 
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4.2 Single Crystal Airfoil Material Constitutive Tests 

Cyclic constitutive tests are being conducted on uncoated primary single 
crystal material (PWA 1480) to develop appropriate constitutive models for 
turbine airfoil applications. Specimens have been oriented in four 
crystallographic directions: <100>, <110>, <123>, and <111>. The specimen 
shown in Figure 2-7(B) has been used for all cyclic tests. To date all cyclic 
tests have been isothermal push-pull type, completely reversed strain cycles. 
Each test condition uses a controlled constant strain rate. 

Three basic types of constitutive data are being obtained: 

1) stabilized hysteresis loops obtained for several different strain 
ranges; 

2) strain rate sensitivity data; and 

3) creep and relaxation behavior obtained during hold periods at various 
points on a given hysteresis loop. 

A typical strain rate sensitivity test is illustrated in Figure 4-10. Figures 
4-11 and 4-12 show results of a typical test in which creep and relaxation 
data are obtained. Figure 4-11 shows the locations on the hysteresis loop at 
which hold periods are imposed. Figure 4-12 shows the decay in the load or 
strain during each of the hold periods. The "negative" stress relaxation, 
shown in the curves A and E, was commonly observed in tests to date. This 
occurs when the stress level is lower than the back stress. After each hold 
period, the cycling is continued until the hysteresis loop again stabilizes. 

In all tests to date, the loop stab! i zed to its prior shape after only a few 
cycles. 

Figures 4-10 through 4-12 are the analog traces taken during testing. The data 
are also being recorded digitally throughout the testing so that virtually 
every loop and hold period is retrievable for analysis. Appendix E shows a 
sample of the digital data for a typical test. 

A total of 15 cyclic constitutive tests have now been completed. The specimen 
numbers of tests at each temperature and orientation are summarized in Table 
4- IV. 

Table 4-IV 


Tests Completed as of December 31, 1984 
Temperature - °C (°F) 


Orientation 

427 (800) 

760 (1400) 

871 (1600) 

982 (1800) 

<100> 

JA61 

JA44 


JA66 

<110> 

KA27 

KA26 


KA22 

<111> 

LA66 

LA67, LA63 

LA68 

LA64, LA69 

<123> 

MA26 

MA25 


MA23 
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Figure 4-12 Typical Relaxation and Creep Data at Various Points on a 
Hysteresis Loop 


Appendix D presents a detailed list of the loading history for each specimen. 
The test conditions are presented in the order in which they occurred. 

The cyclic stress strain data are being reviewed to provide a general 
characterization of the material behavior which can serve as a guide to model 
development and the remaining tests. The general trends observed to date are 
discussed below. 

Figures 4-13 and 4-14 show typical hysteresis loops for all four crystal 
orientations at 982°C and 427 C (1800°F and 800°F) respectively. The strain 
rate shown here is 10 _ 3 in/in/sec. A strong orientation dependence is 
obvious at both temperatures. In general, the amount of inelasticity is seen 
to increase^for orientations further away from the crystal axes, <100>. At 
982“C (1800' , F) the <110> and <123> specimens have nearly identical hysteresis 
loop shapes. These two orientations also have the same apparent modulus ( i . e . , 
1/S 1 33 where S’ 33 is the component of the compliance matrix relating the 
stress and strain in the tensile direction); therefore, the inelastic portions 
of their response should be nearly identical. Construction of a 
stress-inelastic strain curve for the 982°C (1800“F) test. Figure 4-15, shows 
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this to be true. Figure 4-15 was obtained from a classical cyclic 
stress-strain curve which joins the tips of the stabilized hysteresis loops 
for various strain amplitudes. The cyclic stress strain curves for the 982°C 
( 1800° F ) tests are shown in Figure 4-16. The apparent modulus values used were 
taken from the elastic region o of the o first loading portion of each test. 

Figure 4-14 shows that at 427°C (800°F), the <110> and <123> orientations do 
not display such a close similarity as was observed at 982°C (1800°F). This is 
shown more dramatically in the stress-inelastic straiji plot of Figure 4-17. 


PWA 1480 STABILIZED LOOPS 

1 800°F, 10' 3 IN/IN/SEC 


STRESS. KSI 



< 1 1 1 > 
<110> 

< 100 > 


Figure 4-13 Orientation Dependence at 982°C (1800“F) 
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Stress, ksl 


Figure 


PWA 1480 




Stress-Strain Curve 
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Figure 4-16 Cyclic Stress-Strai n Curve at 982°C (1800°F) Constructed of Tips 
of R e : -1 Hysteresis Loops at 10“ 3 in/in/sec 
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Figure 4-17 Stress Versus Inelastic Strain Curve from Cyclic Stress-Strain 
Curve 


In general it was observed that the orientation dependence diminishes at 
higher temperature. This is consistent with the notion that 
noncrystallographic, thermally activated processes dominate at elevated 
temperatures. Another manifestation of the decreased sensitivity to 
orientation at elevated temperatures was the behavior of the specimens on 
initial loading compared to the stabilized behavior. At 427°C and 760“C (800“F 
and 140(rF) slip bursts were observed during the initial cycles of loading as 
illustrated in Figure 4-18. This effect was not observed at 982“C (1800“F). 

The effect was most apparent in the <110> and <123> specimens. These 
orientations have fewer highly loaded slip systems than do the <100> and <111> 
specimens so that slip is more likely to be concentrated. 
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Figure 4-18 Cyclic Hardening Response 


The slip bursts are observed at both tensile and compressive ends of the loop. 
(This behavior is not unlike the Portevin-LeChatelier effect frequently 
observed in superal loys. ) Upon loading in either tension or compression, a 
critical stress is reached at which the material slips very rapidly. When the 
material slip burst strain rate is faster than the controlled strain rate at 
which the test is being conducted, the applied load is reduced until the 
controlled strain rate is again achieved. This results in the small bumps or 
serrations in the loop as shown in Figure 4-18. In general, the stress at 
which these slip bursts occur increases with each subsequent cycle so that a 
"hardened", stabilized hysteresis loop is achieved after several cycles. The 
"hardening" is manifested as both an increased stress range and a decreased 
inelastic strain range (loop width at zero stress). Specimens showing this 
behavior also show easily observable macroscopic slip bands on the gage 
section. 
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The tests conducted to date included an evaluation of strain rate sensitivity 
with emphasis placed on the higher temperature tests. Figure 4-19 shows the 
982°C (1800°F) strain rate sensitivity as measured by the stress amplitude for 
a given strain range. The data have been normalized to the stress amplitude 
observed at 10" 3 in/in/sec. It is clear that the rate dependence is a 
function of orientation and strain amplitudes, with the effect being more 
pronounced for higher strain amplitudes and deviation of the <100> 
orientations. Similarly, the tension-compression assymetry at 982°C (1800“F) 
was found to be rate dependent and also a function of orientation and strain 
amplitude as shown in Figure 4-20. In the limits of elastic loading there 
should be no orientation or strain amplitude dependence. It should be noted 
that in all the tests to date, the initial loading was compressive. Future 
tests will include load cycles which produce tensile yielding first. 
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Figure 4-20 Tension-Compression Assymmetry, Strain Rate Sensitivity 


4.3 Single Crystal Fatigue Tests 
4.3.1 Test Facility 

The test facility used for isothermal and thermo-mechanical fatigue (TMF) 
tests consists of a servo-controlled, closed loop hydraulic testing machine 
with MTS controllers, a low frequency (10 kHz) 20 kW TOCCO induction heater, 
and an Ircon u model 7000 radiation pyrometer, calibrated over a temperature 
range of 260°C to 1371°C (500"F to 2500°F), for temperature measurement. 
Induction heating was selected to accommodate MTS external extensometry and 
specimen inspection and to provide adequate nonisothermal test heating rates. 
The quartz rods of the MTS extensometer, which define a 2.54 cm (1.0 in) gage 
section, are spring loaded against the specimen and did not show any signs of 
slippage during testing. A typical test setup is illustrated in Figure 4-21. 


4-26 



ORIGINAL PAGE IS 
OF POOR QUALITY 



Figure 4-21 Thermo-Mechanical Fatigue Test Rig 


4.3.2 Strain Measurement 

For this research effort, there are two uniaxial deflection measurement 
devices (Figure 4-22) available for elevated temperature fatigue testing of 
hollow specimens: 1) an internal extensometer developed at Pratt & Whitney in 
the 1960s, and 2) a recently developed MTS external extensometer. As 
illustrated in the figure, the internal extensometer requires machined 
internal ridges located at the extremities of the gage section. 

Experience with coated specimens controlled via the internal extensometer 
indicates that thermo-mechanical fatigue (TMF) crack initiation is insensitive 
to specimen internal geometry since cracking predominantly originates on the 
coated surface (the O.D.). However, there is a concern that cracking may 
originate from the internal ridges of uncoated specimens, and that the 
associated crack initiation life of such tests would be indeterminate. Thus, 
it is advantageous to remove the internal ridges and rely on external 
extensometry for deflection measurement. 
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Figure 4-22 Extensometry 

To substantiate the measurement accuracy of the external extensometer, a 
comparison study is being performed between it and the internal extensometer. 
Deflections from each device are recorded from the same specimen by using the 
"piggy-back" setup of Figure 4-22. Previous experience from a single PWA 1455 
( B1900+H f ) TMF test indicated that significant differences between the two 
extensometers occurred only when large cracks were present. During that stage 
of fatigue life, the external extensometer recorded higher deflection ranges 
than the internal. This was interpreted as a limitation of the external 
extensometer because it measures a strain which can be influenced by the 
circumferential location of the dominant crack. The internal extensometer is 
considered to be insensitive to crack location. For the present study, an 
initial set of six coated PWA 1480, <100> orientation, single crystal 
specimens were run at the conditions shown in Figure 4-23. At each test 
condition, specimens coated with each coating type were tested. Observed 
deflection results of these tests (Table 4- V ) indicate that the external 
extensometer consistently measures less mechanical strain than the internal. 
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Figure 4-23 Initial Thermo-Mechanical Fatigue Test Conditions 
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Table 4-V 


Comparison Summary of Strains Obtained from Internal and 
External Extensometry During Thermo-Mechanical Fatigue Cycling 


Strain Ranges (Percent) 


Spec 

ID 

Test 

Condition 
°C fF) 

Coat 

Type 

Thermal 
Ext Int 

Total 

Ext 

Int 

Mechanical 
Ext Int 

JB-9 

427-1038 

(800-1900) 

Out-Of-Phase 

286 


External Extensometer Quartz Rods 
Improperly Located On Specimen 


JB-19 

427-1038 

(800-1900) 

Out-Of-Phase 

273 

0.96 

0.94 

0.13 

0.105 

0.83 

0.84 

JB-10 

427-1038 

(800-1900) 

Out-Of-Phase 

286 

0.97 

0.96 

0.47 

0.42 

0.50 

0.54 

JB-22 

427-1038 

(800-1900) 

Out-Of-Phase 

273 

0.95 

0.93 

0.43 

0.36 

0.52 

0.57 

JB-11 

427-1038 

(800-1900) 

In-Phase 

286 

0.97 

0.96 

1.76 

1.76 

0.79 

0.80 

JB-29 

427-1038 

(800-1900) 

In-Phase 

273 

0.97 

0.97 

1.72 

1.76 

0.75 

0.79 


To better understand the observed deflection differences, parametric tests 
were conducted on uncoated, strain gaged PWA 1480 <100> and PWA 1455 
(B1900 + Hf) specimens at room temperature. The strain gage locations are shown 
in Figure 4-24. For these tests, it was assumed that the average strain gage 
measurement was representative of the true O.D. surface strain and could be 
used as a reference with which to compare the two extensometers. 

The results summarized in Table 4-VI indicate that the external extensometer 
measured strains more consistent with those of the strain gages. 


4-30 



Axial Location 



Gage 


ORIGINAL PAGE IS 

OF POOR QUALITY 


Circumferential 

Location 

3 



Section A- A 


Figure 4-24 Strain Gage Locations on PWA 1480 and PWA 1455 Specimens 


Table 4-VI 

Comparison Summary of Strains Obtained from Internal and External 
Extensometry and Strain Gages During Room Temperature Cycling 

Strain Ranges (Percent) 


Spec 

Material 

External 

Extensometer 

Location* 

PWA 

1480 

2 

PWA 

1480 

1 

PWA 

1480 

1 

PWA 

1455 

4 

PWA 

1455 

3 

PWA 

1455 

3 


* Refer to Figure 4-24 
** Corrected for calculated 


Internal 

0.324 (0.309)** 
0.324 (0.309) 
0.330 (0.315) 
0.188 (0.177) 
0.188 (0.177) 
0.188 (0.177) 

distortion 


External 

Average of 
Strain Gages 

0.300 

0.298 

0.293 

0.296 

0.300 

0.302 

0.178 

0.178 

0.177 

0.173 

0.178 

0.178 
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From an elastic two-dimensional finite element analysis of the specimen 
geometry shown in Figure 2-7(C), it was concluded that both extensometers 
accurately measure deflections at their respective locations. Further, it was 
concluded that the observed differences in strain (at least at room 
temperature) are caused by internal ridge rolling due to the Poisson effect 
and the specimen thickness increase away from the gage section. When the 
internal extensometer strain is adjusted to account for these effects (refer 
to the numbers in parentheses in Table 4-VI), the difference is less than one 
percent for the PWA 1455 specimen. Since the Poisson ratio of PWA 1480 <100> 
is greater than PWA 1455, greater ridge deflection will result. Thus, the 
larger differences of the PWA 1480 specimen can also be explained. Additional 
finite element analysis using PWA 1480 anisotropic properties will be 
performed to verify this conclusion. 

In light of these findings, there is currently no technical reason why the 
external extensometer cannot be relied upon to measure hollow specimen 
strains; however, the decision to remove the internal ridges from future 
specimens and rely solely on the external extensometer is pending the results 
of additional testing and analysis at elevated temperature conditions. 


4.3.3 Single Crystal Fatigue Tests 

Isothermal and thermo-mechanical fatigue (TMF) tests are being conducted to 
define the crack initiation life of coated PWA 1480 single crystal material 
and to provide baseline data for life prediction model evaluations. To 
preclude specimen geometry effects, all fatigue tests (isothermal and TMF) 
will be conducted on hollow specimens, Figure 2-7(0. Initial testing will be 
limited to key variables considered relevant to the creep- fatigue life 
prediction problem. The variables include crystallographic orientation, strain 
range, strain ratio (R e ), strain rate, strain dwell periods, and 
temperature. Both coated and uncoated specimen tests are planned to determine 
whether uncoated testing is useful for evaluating coated crack initiation 
life. The applied coatings are of two generic types which are currently used 
for airfoil oxidations protection: 1) an aluminide diffusion coating (PWA 
273), and 2) a NiCoCrAlY overlay coating (PWA 286). In an effort to 
characterize the material deformation mechanisms which are active during 
fatigue, several specimens will be selected for transmission electron 
microscopy (TEM) examination. 

To date, fully-reversed (R ? = -1) TMF tests of coated PWA 1480 <100> have 
been conducted at the conditions shown in Figure 4-23. These data are 
summarized in Table 4-VII. Cracking information was obtained by acetate 
replication of the specimen gage section at various cyclic intervals. Whenever 
possible, the crack which ultimately led to specimen fracture was identified 
and traced to its initiation site. 
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Separation life 



For the out-of-phase TMF tests, the two coatings were observed to initiate 
cracks differently. The diffusion coating tended to generate circumferential 
cracks, whereas the overlay initiated a multitude of small thumb-nail like 
cracks. Representative fracture surfaces of these specimens are presented in 
Figures 4-25 and 4 - 26 . The concentric cracking of the diffusion coating is not 
unexpected since its fracture ductility is low at temperatures below 649 J C 
(1200 J F). From the replica information (Figure 4-27) it is apparent that lower 
test strain ranges are required to obtain relevant cracking data for this 
coati ng type . 

In-phase TMF crack initiation of coated specimens is less distinct. The 
tensile portion of the test occurs at high temperature; therefore, the 
coatings tend to rumple which makes replication difficult. As can best be 
determined, in-phase testing generates many small cracks early in the fatigue 
life of the specimen; however, these cracks do not grow as rapidly as those of 
out-of-phase testing (Figure 4-28). As a result, the fatigue lives of in-phase 
TMF tests are significantly longer than those of out-of-phase tests. 

A life correlation of the fatigue tests is presented in Figure 4-29. For 
out-of-phase tests, the fatigue life of the diffusion coated specimens is 
shown to be more sensitive to total strain range than that of the overlay 
coated specimens. Thus, at strain ranges below about 0.5 percent, the 
diffusion coated specimen lives are expected to be higher than those of 
overlay specimens. 

A general matrix of PWA 1480 fatigue test conditions planned for this task is 
provided in Table 4-VIII. 
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(A) CIRCUMFERENTIAL CRACKING 


~6X 


50X 


B) FAINT FATIGUE STRIATIONS 



Figure 4-25 PWA 273 Coated PWA 1480 <100>; Specimen JB-22; After Being 
Thermo-Mechanical ly Fatigue Tested at 427°C - 982“C (800°F - 
190Q J F) Ae-t = ^0.275 Percent, 1 CPM, Out-Of-Phase for 3772 
Cycles 
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~10X 

A) MULTIPLE SMALL THUMB-NAIL LIKE FATIGUE CRACKS 



40X 

B) FAINT FATIGUE STRIATIONS 


Figure 4-26 PWA 286 Coated PWA 1480 < 100 >; Specimen JB-9; After Being 

Thermo-Mechanical ly Fatigue Tested at 427 ^C - 982 J C (800 F - 
1900°F ) Aet = ^0.4 Percent, 1 CPM, Out-Of-Phase for 1878 Cycles 
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Figure 4-27 Coated PWA 1480 <100> Out-Of-Phase Thermo-Mechanical Fatigue 
Replica Data 



~10X 

Figure 4-28 PWA 273 Coated PWA 1480 <100>; Specimen JB-29; After Being 
Thermo-Mechanical ly Fatigue Tested at 427"C - 982°C (SOO^F - 
1900°F) Aet = +0.4 Percent, 1 CPM, In-Phase for 10340 Cycles. 
Coating rumpling and many small cracks are the dominant 0.0. 
surface features. 
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Thermo-Mechanical 
PWA 1480, 427°C - 


Fatigue Life Versus Strain Range for Coated 
1038 C (800“F - 1900 d F) 


4-38 



Table 4-VIII 


Summary of Planned Task III PWA 1480 Fatigue Tests 



Notes: 1: 1 - <100>; 2 - <110>; 3 - ♦ - <123> 

2: - - uncoated: Cl - PWA 286 NICoCrAlY overlay; C2 ■ PHA 273 aluelnlde diffusion 
3: Cephas Is will be on conducting the coated specimen tests 
4: Special Cycle 
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SECTION 5.0 


SUMMARY OF RESULTS 


During the first year of effort, the following tasks were accomplished and 
preliminary results obtained. 

o The materials selected for the program are as follows: single crystal 
alloys, PWA 1480 and model Alloy 185, an overlay coating, PWA 286, 
and aluminide coating, PWA 273. All but the Alloy 185 were selected 
because they are typical of the most advanced materials being used in 
gas turbine blades. Alloy 185 was selected because its properties 
differ greatly from those of PWA 1480 and will be used to test the 
range of applicability of the constitutive and life prediction models 
developed in this program. 

o Single crystal specimen orientations were chosen to be <100>, <110>, 
<111> and <123> to represent the full range of slip systems and slip 
system interaction. 

o Physical and thermal property tests on PWA 1480 were completed at 
Southern Research Institute. 

o Specimen preparation is well underway. Casting of all PWA 1480 
specimens is complete, and casting of Alloy 185 is underway. 

Machining of these specimens is in process. Coating specimen 
fabrication is being paced by the success or failure of exploratory 
test procedures. 

o Literature searches were completed on related constitutive and life 
prediction models, and are included in the appendices. 

o The Walker unified viscoplastic constitutive model has been extended 
to all the slip systems of a face centered cubic crystal material. 

The initial mathematical development has been completed and 
programmed and is being used to correlate initial laboratory data. 

o Walker's isotropic constitutive model has been successfully used to 
correlate limited data obtained for overlay coating, PWA 286. Model 
constants derived compare reasonably with those from other high 
temperature alloys. 

o Thermal property tests were completed for PWA 1480. Properties 

included conductivity, coefficient of thermal expansion and specific 
heat. 

o Monotonic tensile and creep data were obtained for PWA 1480 and 
overlay coating PWA 286. 
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PWA 1480 constitutive data were obtained on fifteen of forty 
specimens. The constitutive models and their constants are being 
evaluated with these data. From this study the test conditions for 
the next round of constitutive tests will be determined. 

o A modern external MTS extensometer has been evaluated relative to the 
older and well-established Pratt & Whitney internal extensometer and 
was found to function accurately and reliably in thermo-mechanical 
fatigue (TMF) testing. Use of the external extensometer in this 
program will eliminate the requirement for internal ridges on the 
test specimens. 

o Isothermal and TMF testing of coated and uncoated PWA 1480 has 

started. Initial out-of-phase cycle TMF tests indicate that at strain 
ranges greater than 0.5 percent, the overlay coated specimens have 
longer lives than do those coated with an aluminide. The reverse is 
true below the 0.5 percent strain range. 
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SECTION 6.0 


FUTURE WORK 


During the second year of the base program, the following work will be 
accomplished: 

o Specimen preparation will be completed, including the machining of 
PWA 1480 specimens, casting and machining of the Alloy 185 specimen, 
and complete fabrication of all coating specimens. 

o Physical property tests will be run on the two coating materials, PWA 
286 and PWA 273. 

o Level I PWA 1480 experiments will be completed including forty 
constitutive tests and forty fatigue tests. 

o Development of the constitutive and life prediction models will 

continue in parallel with the above tests. Two models of each will be 
selected for further evaluation in Task V. 
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APPENDIX A 

EXPERIMENTAL DATA GENERATED FOR 
COATING CONSTITUTIVE MODEL EVALUATIONS 
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Figure A-l 533°C (100Q°F) PWA 286 Bulk (Hot Isostatically Pressed) 
Relaxation; Unexposed 
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Figure A-2 760°C (1400°F) PWA 285 Bulk (Hot Isostatically Pressed) 

Relaxation; Unexposed 


A-3 


Stress, psi Strain, in/in 




Time, minutes 

Figure A-3 871°C (1600°F) PWA 286 Bulk (Hot Isostatically Pressed) 

Relaxation; Unexposed 
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Figure A- 4 982°C (1800°F) PWA 286 Bulk (Hot Isostatically Pressed) 

Relaxation; Unexposed 
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APPENDIX B 


LITERATURE SURVEY - SINGLE CRYSTAL AIRFOIL MATERIALS CONSTITUTIVE MODELS 


This appendix contains a review of constitutive models or modeling approaches 
which may be applicable to single crystal and transversely isotropic 
materials. The discussion includes work which has been reported in open 
literature and likely extensions of this work. 

Macroscopic Continuum Theory Formulations 

In 1948, Hill (Reference 1) proposed a generalization of the von Mises yield 
function for states of anisotropy possessing three mutually orthogonal planes 
of symmetry. Hill's theory remains the most widely used formulation and is 
currently coded into general purpose nonlinear finite element programs such as 
MARC and ABAQUS. 

Various modifications to Hill's original theory were published in the 
literature between 1948 and 1977. Theories presented by Edelman and Drucker 
(Reference 2), Yoshimura (Reference 3), Williams and Svensson (Reference 4), 
and Baltov and Sawczuk (Reference 5) are basically modifications to Hill's 
theory to account for Bauschinger effects associated with kinematic hardening. 
Further work on anisotropic plasticity may be found in References 6 through 
11. Hill's original theory allowed the anisotropic yield surface only to 
expand, whereas the later modifications allowed the anisotropic yield surface 
to both expand, translate, and change shape in stress space. However, the 
modifications to Hill's theory employ kinematic shift tensors, a-jj , which 
depend on plastic strain. The various shift tensors therefore vanish in the 
initial state of the material in which the plastic strain is zero and do not 
reflect a difference between tensile and compressive loading as observed in 
experiments with nickel-base single crystal materials. In effect the initial 
state is governed by Hill's theory. 

The yield surface employed by Edelman and Drucker may be written in the form: 

f = \ M ijkl (s ij " n,c ij ) (s kl “ mc kl ^ ^ 

where Si j is the deviatoric stress tensor, 
c-jj is the plastic strain tensor, 
k 2 denotes the size of the yield surface, and 
M ijkl ls a fourth order anisotropy tensor. 
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If f < 0, the material is elastic. If f = 0, the yield condition is reached 
and the stresses are on the yield surface. In this theory, the kinematic shift 
tensor, ajj, has the form a-M = mcij and vanishes in the initial virgin 
state of tne material. This tensor translates the yield surface in stress 
space and therefore introduces kinematic hardening to account for the 
Bauschinger effect. 

Yoshimura employed the yield function: 

f ■ 1 M ijki s u s ki - Bs ij c u - k2 - 


in which the yield surface translates through the term BsjjCjj, which also 
vanishes in the initial virgin state of the material, while that of Williams 
and Svensson has the form: 


f - (I ijki +Ac ij c ki> fij'S-j (L 


pqrs a pq CT rs 


->}{ 


a, i + Ci i (L a a 

kl kl pqrs pq rs 


-m) j- - 


The yield function employed by Williams and Svensson was generalized in 1977 
by Lee and Zaverl (References 12 and 13) whose yield surface may be written in 
the form: 


f 


I M ijkl 



“ij* (a kl “ “kl* 



This form of the yield surface will now be considered in more detail since 
this type of formulation is readily adapted for use with the unified 
viscoplastic models whose employment will be required in the description of 
the time-dependent anisotropic material behavior of single crystals and 
directionally solidified materials at elevated temperature. Several aspects of 
a time- independent formulation will be discussed, including 
tension/compression asymmetry, flow and hardening rules, the formulation of 
the elastoplastic stiffness matrix, and finally, the notion of 
incompressibility. Before proceeding, it will be convenient to adopt a Voigt 
rotation for the tensor indices so that the fourth rank tensor can be 

written in matrix form. To this end, we introduce the contracted notation for 
the stress components: 


r l =CT ll ! 


ct 2 = ct 22 1 




3“ 33 5 


CT 4 =cr 12 =CT 21 ' 0 5 =0 23 =u 32’ °6~ u 13" u 31’ 


CSr-0 


a c =Oi 0 =a. 


and the plastic strain components are: 


c r c ll* c 2 =c 22’ 


c 3 =c 33 * 


c 4= 2c 12= 2c 2i> 


c 5 =2c 23 =2c 32 s 


C 6 =2c 13 =2c 31 ' 


Any component of the tensor Mijkl ma Y then be written as a two index matrix 
component. For example, M 12 32 = M 45 an< * ^11 22 = ^12» etc. 

T T TT 
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Lee and Zaverl's yield surface may now be written in Voigt index notation as: 


f = H 13 (a, - a,) U } - aj> - k 2 . (B-l) 

If the material has three orthogonal planes of symmetry (such as the three 
crystallographic faces in a face centered cubic nickel -base superalloy) and 
the principal axes of anisotropy coincide with the Cartesian reference axes, 
the matrix M-f j can be written from symmetry arguments in the form: 


Mil 

Ml 2 

m 13 

0 

0 

0 


M 21 

M 2 2 

M 2 3 

0 

0 

0 


M31 

M 32 

^33 

0 

0 

0 

(B-2) 

0 

0 

0 

M44 

0 

0 


0 

0 

0 

0 

^55 

0 


0 

0 

0 

0 

0 

M66 



Equation (B-l) may also be written in the more convenient form: 

f = - L ^ a - - J, ( B-3 ) 


where 


K. = 2M..a. and J = K 2 
i U J 


M . .a.a . . 
1J 1 J 


(B-4) 


Again, by symmetry, the vector Lj may be written in component form as: 

CL] = ll v L 2 , L 3 , 0, 0, 0]. (B-5) 


A physical interpretation of Mjj, Lj and J may be made in the following 
manner. Let Y+i, Y+2, Y+3 denote the tensile yield stresses in the xi, 

X2, X3 directions and let the absolute values of the compressive yield 
stresses in these directions be denoted by Y_i, Y_2, Y_3, respectively. 

From equation (B-3), with f=0, we may write for simple uniaxial tension and 
compression tests in the xi direction: 

M 11 Y +1 " L 1 Y+ 1 = J in tens1on > (B-6) 


and 


M n v ‘i * 4 Y -i 


in compression. 


(B-7) 
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(B-8) 


Equations (B-6) and (B-7) show that: 

M 11 s Y^Y“ 5 M 22 = Y + 2 Y _2 * M;33 = 


We are at liberty to choose J=1 in equation (B-3) since this merely scales the 
values associated with M-jj and L-j. Hence, 


'll * 'V-l' -1 • M 22 ■ <V-2 rl • M 33 ' (V-3 1 ' 1 • 


(B-9) 


and the components of M-jj are determined from measurements of the tensile 
and compressive yield stresses along the three crystallographic directions. In 
pure shear equation (B-3) may also be shown to give: 


m 44 - <y 4 ) 


-2 


M 55 = ( V 


-2 


'66 


( V 


(B-10) 


where Y4, Y5 and Yg are the torsional yield stress associated with the 
stress components 012, 023 and 0^3, respectively. By substituting 
equation (B-9) into equations (B-6) and (B-7), we find that: 



U = 



( B— 1 1 ) 


which determines the values of the vector in terms of the measured yield 
values in tension and compressor. Thus, initial tensile and compressive yield 
stresses along the crystallographic axes, in addition to torsion tests in the 
12, 23 and 31 directions, determine the initial values of the matrix M-jj and 
the vector Lj. 

To complete the rate- independent formulation it is necessary to specify a 
suitable flow rule and to choose suitable hardening rules which determine the 
subsequent hardening evolution of the variables 04 and k 2 in the 
specification of the yield surface in equation (B-3). the equivalent stress, 
k, is defined by the yield condition in equation (B-3), viz: 

k = CM pq (a p • V {a q ■ “q 3 ^ * (B ' 12) 

and the equivalent plastic strain increment, dR, may be defined through the 
work relation: 


kdR = a. dc.. . 


(B-13) 
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If the yield function and the plastic potential are assumed to be identical, 
the flow rule is obtained by the usual relationship: 


dc. = x~- 
1 3 a.. 


( B-14) 


The non-negative function \ is obtained from the usual consistency relation 
such that when the stress state is on the yield surface. 


df = da. + Jp- dM. . + da . + dk 

3a i 1 3 ^ij 3a i 1 3 * 


0 . 
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The constitutive relation may be written in the form: 


do, = 0,. d tj , 

where the elastoplastic matrix D*j is given by: 

k n 3 f 3 f 


* i 

D. . = D. • - 4 r 

10 10 a a_ 3T 


m 


da. 


D. -LL — — D . , 

ip 3 a p 3 a q qo » 


m 
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(B-17) 


wi th 


2k iJl + p- 

sTT a_ 3 f 
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3 a, 


lf_ D .. 1L 

3a. u 10 3a- 

! J 


3f 3M 1j 3f 3ai 

TFT7. TR 3a. TR - 

* J ' 
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In this relation for d, it is assumed that Mjj, aj and k are functions of 
the cumulative plastic strain, R, and the matrix in equation (B-17) is 
the anisotropic elasticity matrix: 


Dll 

d 12 

d 13 

0 

0 

0 

°21 

d 22 

d 23 

0 

0 

0 

d 31 

d 32 

d 33 

0 

0 

0 

0 

0 

0 

d 44 

0 

0 

0 

0 

0 

0 

°55 

0 

0 

0 

0 

0 

0 

d 66 


( B- 1 9 ) 


when referred to the crystallographic symmetry axes. 
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We now examine the consequences of assuming that the plastic strain is 
incompressible. Setting dcj + deg + dc3 = 0 in equation (B-14) yields 
the relations: 


Mu + ^12 + M 13 = 0 , (B-20 

^21 + ^22 + ^23 = 0 , (B-21 

M31 + M32 + M33 = 0 , • ( B- 22 

Ll + L2 + L 3 = 0 . (B-23 

In a face centered cubic material, such as a nickel-base superalloy single 
crystal, the properties along each of the crystallographic axes are identical, 
so that. Mu = M22 = M33 and M12 = M23 = M13. 

From the preceding discussion, it is evident that the anisotropy matrix, when 
specialized to the cubic symmetry appropriate for single crystals, may be 
written as: 


In t"h 


~¥ll 


’ll T"ll 


[M] = -^n tMji 


while the elasticity matrix has the form: 


Dll 

Dl2 

°12 

0 

0 

0 

°12 
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d 12 

0 

0 

0 

D i2 

d 12 

D 11 
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0 
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D 44 
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0 
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In order to more clearly investigate the ability of this type of theory to 
reproduce the orientation dependence of the initial yield stress, we will 
assume that the material behaves as an anisotropic elastic-perfectly plastic 
material. Consequently, we may set a-j = 0 for the present discussion, since 
the kinematic shift vector aj changes only with plastic deformation and we 
are concerned only with the prediction of the tensile yield stress. Also, 
since we are interested only in the initial yield stress and not on the 
subsequent tensile hardening of the stress-strain curve, we may set: 

M-jj = constant; k = constant; and ai = 0. 

Consider a uniaxial tensile bar specimen of the single crystal material 
oriented along the direction ON in Figure B-l. When reduced to one-dimensional 
form, the elastic-perfectly plastic constitutive relation of Lee and Zaverl 
appropriate to tensile loading reduces to the following relations: 

a= E(e,*) for e < Y + (e,0)/E(e,«) , (B-26) 

a = Y + (e,tf) for e > Y + (e,*)/E(e,0) . (B-27) 


In these expressions E(e,^) is the value of Young's modulus in the direction 
ON given by the expression: 


E(M) * 

with 


D ll +D 12 . / 1 2 \ / 2 2 . 22 . 22, 

{ °1 i + 2Di 2 > ( i- D i2 > \ D 44 D H“ D 12 ) xy xz yz 


, (B-28) 


n x = -sin 0; n y = sine cos$; n z = cose cosip . (B-29) 

The tensile yield value in the direction ON is given by the relation: 


Y + (e,W = 


' Y + -Y_ 


fY + -Y 

VC 


4 {y 7T + x(n x n y + n x n z + n y n M 


2 {yJ-_+ ^ n x n y + n x n z + n y n M 


, (30) 


with 


\ = 4 M 44 + 2Mi2 - 2Mn = 4 M 44 - 3M n . 


(B-31 ) 
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The constant A measures the deviation from isotropy and 0 in the case of 
an isotropic material, while the constants Y+ and Y_ are the tensile and 
compressive yield stresses along any of the three crystallographic axes, the 
properties along each axis being assumed to be equivalent. 


2 [ 001 ] 



Point E In the Stereographic 
Triangle ORS Represents the 
/ yX Direction ‘ON 1 In the Cubic 
\ Lattice Shown Above 

• E \ 

_/u R] 

[001] — [Oil] 

(0 * 0 °, * « 0 °) (« ■ 45 °, 0 » 0 °) 


Figure B-l A Single Crystal Bar Oriented Along ON at Angles e and $ With 

Respect to the Crystal Axes X, Y, Z is Located at Point E in the 
Stereographic Triangle ORS. 


Figures B-2 and B-3 show actual and predicted variation of the tensile yield 
stress, Y+(g,^), as a function of the orientation angles e and $ for the 
case of single crystal PWA 1480 at 593°C (1100°F). It is evident from Figure 
B-2 that the theoretical prediction of the yield stress variation with e for 
the case where i M°, i.e., along the [001]-[011] boundary of the stereographic 
triangle in Figure B-3, shows a yield stress which decreases continually from 
a maximum at e=0° (along the [001] direction) to a minimum at e=45° (along the 
[Oil] direction). A further decrease in the tensile yield stress occurs as 4* 
increases towards the [001]-[111] boundary of the stereographic triangle. An 
absolute minimum in the o theoretical prediction occurs in the [111] direction 
where 9=45° and 35.26°. In the experimental data shown in Figure B-3, the 
variation of yield stress along the [001]-[011] boundary of the stereographic 
triangle shows a minimum value around e=35°. 
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[ooi] 8“*" [on] 

Figure B-2 593°C (nOO°F) Tensile Iso - Yield Contours for PWA 1480 Based 

on Lee and Zaverl's Theory; MPa (ksi) 



Figure B-3 Experimental Tensile Iso - Yield Contours for PWA 1480 





The experimental data are more consistent with the notion that the material 
yields when the resolved shear stress reaches a critical value in one of the 
slip directions on one of the crystal! ographic slip planes. For slip which 
occurs in the jlllf octahedral planes in the <110> type slip directions, the 
variation of yield stress with angle e along the [001]-[011] boundary of the 
stereographic triangle is given by the relation, 


Y + (M=0°) 


Y + (e=0\ $= 0 °) 
cos q (sine + cose) 


(B-32) 


This expression is plotted in Figure B-4 and exhibits a "U" shaped curve in 
which the minimum yield stress occurs at e = 22.5°. The experimental data are 
also plotted in this figure for comparison purposes. More will be said 
concerning the micromechanics of crystalline slip in the next section. The 
point to be made here is that a minimum is observed in the experimental data 
which is more consistent with a crystalline slip theory than with a continuum 
formulation. 



Figure B-4 The Experimental 593°C (1100°F) Yield Strength of PWA 1480 Vs. 
Orientation Along the [001] - [Oil] Boundary of the 
Stereographic Triangle 
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In general, the as-cast single crystal does not exhibit complete symmetry 
along the three crystallographic axes. For example, Kear and Piercey 
(Reference 14) show that creep rupture lives along the [001], [100] and [010] 
crystallographic axes for single crystals of Mar-M200 are 325, 71 and 82 hours 
respectively at a stressjevel of 670 MPa (100 ksi). After heat treatment for 
one hour at 1232°C (2250°F) and 32 hours at 871°C (1600°F), the creep rupture 
lives were 459, 365, and 443 hours, respectively. The heat treatment produces 
equality of behavior along the crystallographic axes. In the as-cast material 
prior to heat treatment, the inequality of material behavior along the 
crystallographic axes is attributed to dendrite morphology. Initially, 
material bars cut parallel to the continuous dendrite arms along the [001] 
growth axes exhibited a longer creep rupture life than material tested 
parallel to the discontinuous side arms of the dendrites along the [100] and 
[010] directions. The significance of this observation is that better 
correlation with the experimental results may be obtained if the cubic 
symmetry restriction is relaxed and the matrix M-j j is assumed to have the 
general form depicted in equation (B-2). 

Directionally solidified materials consisting of columnar grains of single 
crystals oriented at random about the grain growth direction may be modeled as 
transversely isotropic material if a sufficient number of grains are sampled 
in the transverse plane. In modeling single crystals with cubic symmetry only 
two independent constants are needed in the matrix Mjj and three in the 
elasticity matrix D-m. The incompressibility of plastic flow reduces the 
number of independent constants in M-j s from three to two. Five constants are 
required in the matrices D-j j and Mj j tor directionally solidified 
materials which exhibit transversely isotropic behavior. However, the 
incompressibility condition for the plastic strain reduces the number of 
independent constants in Mjj from five to three. These transversely 
isotropic matrices have the form: 

0 12 Dig 0 0 0 

D 11 °13 0 0 0 

D 13 D 33 0 0 0 , (B-33) 

o o d 44 0 0 

o o o d 44 0 

oooo |(d u -d 12 ) 
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and 



M 11 

M 12 

M 13 

0 

0 
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M 12 

M 11 

M 13 

0 

0 

0 

[M] = 

M 13 

M 13 

f1 33 

0 

0 
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0 

0 

0 

M 44 

0 
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0 

0 

0 

0 

M 44 
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0 

0 

0 

0 

0 

i (M ir 

in which X3 is the columnar growth direction of the single cr 
The incompressibility assumption gives: 

M 12 = ‘ 

7 (2M 11 

‘ M 33 ) 

and 

= ■ \ M 33 

5 


which reduces the number 

of independent constants 

in Mj 

j from 


(B-34) 


(B-35) 


The requirement of transverse isotropy can be met in practice if the number of 
grains is large with respect to the reference section (or equivalently the 
grain size is very small), such as the wall thickness of an airfoil. It is 
also necessary that the growth axes of the columnar crystals coincide with the 
growth axes of the composite directionally solidified material and that the 
secondary crystal orientations about the growth axes are distributed at 
random. Usually, limitations imposed by both the solidification process and 
the specimen size make it extremely difficult to meet the preceding 
theoretical requirements and significant scatter must be expected. Due to this 
scatter, it may again be necessary to revert to the more general matrices 
given in equations (B-2) and ( B- 19 ) to fit the actual material behavior. 


A constitutive equation for modeling secondary creep in directionally 
solidified eutectic alloys, based on a generalization of the Bailey-Norton 
creep law to transversely isotropic materials, was proposed by Johnson 
(Reference 16) in 1977. The theoretical model was developed by assuming that 
the constitutive equation is: 


U 


3 <* 

3a ij 


(B-36) 


where i is a function of the three invariants of the deviatoric stress tensor 
s-jj in addition to the invariants I* j s j -f and s i j I j k s ki » 
where ifj = d-jdj with d-j denoting a unit vector along the growth 
direction of the transversely isotropic material. A similar theory, in which 
the Bailey-Norton creep law is replaced with Robinson's constitutive model, 
was proposed by Robinson in Reference 17. The matrix Mjj which relates the 
inelastic strain rate to the stress components (a-jj - fi-jj or a-j,-) is of 
the same form in the models of Robinson, Johnson and Lee and Zaverl , and 
contains three independent constants for the case of transverse material 
isotropy and two in the case of cubic symmetry. 
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Stouffer and Bodner (Reference 18) have extended the Prandtl-Reuss flow rule 
to allow for anisotropy in which the inelastic strain rate is related to the 
deviatoric stress through the relationship: 



where the matrix X-t j is a function of stress and deformation history and is 
constrained to satisfy the incompressibility condition. Further work on this 
model is being pursued by Stouffer. 


Micromechanical 


Theory Formulations 


An analysis of single crystal superalloys undergoing steady state creep 
deformation was presented by Paslay, Wells and Leverant (Reference 19) in 1970 
using a theoretical formulation based on crystallographic slip theory of face 
centered cubic materials. In 1971, the theory was applied by Paslay, Wells, 
Leverant and Burck (Reference 20) to describe the creep behavior of single 
crystal nickel-base superalloy tubes under biaxial tension loading conditions. 
Steady state creep formulations suitable for the analysis of single crystals 
were used by Brown (Reference 21) in 1970 and by Hutchinson (Reference 22) in 
1976 to predict the behavior of polycrystalline materials whose aggregate 
consists of randomly oriented single crystal grains. Recently, Weng (Reference 
23) has developed a single crystal creep formulation which accounts for 
transient (primary) as well as steady state (secondary) creep. However, in 
order to describe the combined plastic and creep behavior of polycrystalline 
materials, Weng combines the rate- independent and rate-dependent components of 
crystal behavior in such a way that each component is governed by a separate 
constitutive relation, i.e., plasticity and creep are assumed to be uncoupled 
phenomena. 


A special issue of the ASME Journal of Engineering Materials and Technology 
published in October, 1984, is devoted entirely to the micro and macro aspects 
of constitutive equations. This edition contains the latest version of a 
physical theory of plasticity and creep proposed by Lin (Reference 24) which 
averages a single crystal model to obtain a stress-strain relationship for an 
isotropic polycrystalline aggregate. Krempl and Lu (Reference 25) also discuss 
the latent hardening effects observed in fully annealed AISI type 304 
stainless steel under biaxial in-phase and out-of-phase strain cycling at room 
temperature. 


In the decade of the seventies, the creep and plastic responses of materials 
were combined into unified viscoplastic formulations (Reference 15). These 
formulations differ from steady state creep theories by introducing history 
dependent state variables to account for primary creep and plasticity. Most of 
these unified theories exhibit an elastic response under instantaneous 
deformation and can be integrated easily using an explicit forward difference 
method. An anisotropic slip theory formulation suitable for representing 
viscoplastic material behavior of single crystal materials at elevated 
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temperature can therefore be obtained by replacing the secondary creep model 
in Paslay, Wells, Leverant and Burck's formulations with a suitable unified 
viscoplastic formulation. The single crystal slip model has been embedded in a 
self-consistent formulation to model the behavior of transversely isotropic 
directionally solified materials. This work is proceeding under a grant to the 
University of Connecticut (Reference 26) sponsored by the NASA-Lewis Research 
Center. Similar work is being undertaken by Stouffer at the University of 
Cincinatti also under sponsorship from the NASA-Lew,is Research Center. 

The constitutive formulation represented in Reference 6 has treated the single 
crystal as if it were a homogeneous single crystal material. In fact, 
nickel -base superalloys such as PWA 1480 are composed of face centered cubic 
y 1 particles of cuboidal shape embedded in a face centered cubic y solid 
solution matrix. Shah and Duhl (Reference 27) have studied the behavior of PWA 
1480 in detail and have concluded that the deformation is strongly influenced 
by the y-y 1 structure. Constitutive equations for the composite y-y 1 mixture 
could be set up in the manner employed by Nemat-Nasser and his colleagues in 
References 28 through 30. The work reported in Reference 30 was sponsored by 
the NASA-Lewis Research Center and can easily be adapted to the modeling of 
y-y' mixtures of single crystal materials. Further references pertinent to the 
development of viscoplastic slip constitutive models can be found in Reference 
26. 
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APPENDIX C 


LITERATURE SURVEY - SINGLE CRYSTAL LIFE PREDICTION 


C.l Introduction 

All of the life prediction models published in literature have been developed 
for conventional multi grained materials which are considered isotropic for 
analysis purposes. It is generally accepted that structural failures result 
from subjecting materials to inelastic strains. Applications to anisotropic 
single crystal material are documented only for a few of these models. The 
role of coatings on cracking is not modeled although coated specimen data are 
used by some investigators. 

Tne published life prediction models relate time to failure (life) to the 
inelastic strain or its time dependent (creep) and time independent (plastic) 
strain components. Three different approaches are used: 

1) life related to macroscopic inelastic strain and cycle parameters 
(phenomenological models); 

2) life related to damage explicitly caused by creep and plasticity 
(cumulative damage models); and 

3) life related to local inelastic strain at a void and/or crack tip 
(crack growth models). 

A number of the existing models which are surveyed in Sections C.2 through 
C.4.2 of tnis appendix nave potential for extension to hot section anisotropic 
materials under the current oase program. However, methods must be developed 
to account for the important effects of the base alloy anisotropy and the 
coatings. In essence, coatings are thin layers of isotropic materials on the 
surface of the base alloy. A review of published metallograpiiic work on single 
crystal crack initiation and the role of coatings on crack initiation is 
included in Section C.5. 

C.2 Phenomenological Models 

The Coffin-Manson model (References 1 and 2) has been used as the starting 
point by many investigators interested in predicting cracking life of 
structures subjected to cyclic loadings. Originally defined for low and 
intermediate temperature application, this model assumed a logarithmic 
relationship between plastic strain range and low cycle fatigue life. 

ae p N f B = C (C-l) 


where Aep = plastic strain range, 

Mf = number of cycles to crack, and 
B, C = material constants. 
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At high temperature conditions, damage due to both plastic and creep strains 
can occur. To account for the creep effects, the Coffin-Manson equation has 
been modified by substituting inelastic strain range ( a c-f ne -j ) for Aep 
(Reference 3). This allows the effects of creep to be implicitly included in 
the model, but only for those load conditions used in material property tests 
for model development. Using Aej ne ] as the single correlation variable does 
not allow extending this model to other load conditions, i.e., cycle shapes, 
hold times, strain rates, mean stresses, and temperature cycles which cause 
different amounts of creep. 

The Frequency Modified Life model of Coffin extends the Coffin-Manson model at 
elevated temperatures by introducing a frequency correction (References 4 and 
5). This model is described by the following equation: 

Ae inel ' {H f ' vk ’ 1)/? * C 2 (C ' 2) 


where v = cycle frequency, and 

k, 13, C2 = material constants. 


Frequency correction is important at high temperatures because of strain rate, 
dwell time at a given stress, and possible environmental interaction effects 
on low-cycle fatigue (Reference 6). Material property data obtained by 
specimen testing at elevated temperatures inherently reflect load cycle 
frequencies used in conducting the tests. The frequency modification provides 
a single correction to extend data to other frequencies than those used in 
testing. However, corrections for load and temperature cycle shape variations 
are not possible by tnis method. 

It is important to note that the frequency correction used in this model 
distinctly depends on temperature. As stated above, a correction is important 
at high temperatures while no correction appears necessary at low temperatures 
where creep cannot occur. This implies that the frequency correction depends 
on the temperature level, analogous to creep rate, and a single low-cycle 
fatigue correlation does not apply at all "high" temperatures (Reference 5). 
However, satisfactory correlations can be achieved for individual isothermal 
conditions (References 5, 7 and 8). 


The Frequency Separation method (References 9 and 10) further improves the two 
previous models by accounting for the unbalanced tensile and compressive 
portions of the load cycle and recognizing that compressive and tensile creep 
damages are different. Tnis model can be described by an equation very similar 
to that of the Frequency Modified model, except the absolute frequency is 
replaced by the tension going frequency term and a frequency unbalance 
correction strain range is used: 



(C-3) 
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(C-4) 


where 4 e] nel 



s 


vj = tensile going frequency (1/tii.ie for e-; ne ^ >0), 

v c = compression going frequency (1/time for'ej ne i < 0), and 

k, ki, p, Cg = material constants. 

This formulation is oased on the concept that low cycle fatigue damage is 
similar to crack growth process where damage occurs mainly when the crack is 
opening as intense local tensile strains are formed and the crack tip advances 
by local fracture. 

A different frequency separation approach is described in Reference 11, using 
a modified equation of the following form: 



Both unbalance corrections account for the relative times in compression and 
tension but again do not consider stress levels for each. For this reason, 
this model is sometimes interpreted to be more appropriate for environmental 
damage than creep effects on low cycle fatigue (Reference 11). 


The Tensile Hysteretic Energy model was developed by Ostergren (Reference 12) 
to account for the observed influence of strain cycle shapes on low cycle 
cracking life. The model relates cracking life to inelastic tensile strain 
energy and in this manner includes the combined effects of both unbalanced 
time and unbalanced loading. This result is yet another formulation very 
similar to that of Frequency Modified Life, but includes the effects of time, 
temperatures and loads. 


where Alff 


aW t 





tensile hysteretic energy. 


(C-5 ) 


1 

t 0 +t T- t C for t T > T C 

v = 

1 

T 0 for T j £ T c 

t q = time per cycle for continuous cycling, 

r T = tension hold time, and 

Tq = compression hold time. 
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The basis for this model also is that low cycle fatigue is considered to be 
essentially a problem of crack propagation from inherent defects or 
microcracks which initiate in a very small number of cycles. Consequently, it 
is assumed that the deformation which contributes to damage by propagating the 
crack occurs during the portion of load cycle when the crack is open. Next, it 
is recognized that both tensile and compressive hold times increase time 
dependent inelastic strains. But tensile hold leads to compressive mean stress 
while compressive hold leads to tensile mean stress and a larger tensile 
portion for the hysteresis loop. Noting that the compressive hold times were 
found to be more damaging than tensile hold times in high temperature 
materials, it followed that the net tensile hysteretic energy is used as the 
measure of damage. 

For computation convenience, it may be assumed in most cases that: 

aW x = a ' a-r Ae • . 

T T inel 

where a = shape factor for tensile hysteresis loop, and 
CT y = a + -Jr = maximum tensile stress. 


Consequently, the model formulation may be changed to: 

°T ' 4e inel ' (N f ' '’ l ‘~ 1) ‘ 3 * C 'l (C - 6 > 

The frequency modification term definition used in this model is empirically 
selected. Modified definitions have been used when environmental instead of 
time dependent strain effects are modeled (Reference 11). In all cases, the 
frequency modification may depend on temperature, however. 

The tensile hysteretic energy model was recently investigated for application 
to B-1900+Hf material under high temperature fatigue loads (Reference 3). The 
results showed consistent under prediction for some data sets, indicating that 
the total tensile hysteresis area is not equally damaging and/or the damaging 
portion of hysteresis area was not accurately sized because effects of back 
stress were not included. A second difficulty was that of accurately 
determining the area when inelastic strains are small, an important condition 
for most long life turbine blades. Both of these results are related to 
modeling the constitutive behavior of the material and reflect the current 
state of understanding of such behavior. 
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The Tensile Hysteretic Energy model was recently also applied to limit 
correlation of isothermal low cycle fatigue (LCF) and thermomechanical fatigue 
(TMF ) data of PWA 1480 single crystal material (References 13 and 14). In 
order to account for the crystal geometry and loading axis orientations, the 
hysteretic energy term was normalized by the resolved octahedral normal stress 
and the modulus in the loading direction. Excellent consolidation of the LCF 
data was obtained for [001] and [111] specimens in the uncoated condition and 
including test temperatures of 427°C, 760°C and 980°C (800°F, 1400°F and 
1800°F) . 

The TMF testing involved a very limited test of PWA 273 coated specimens of 
[001] and [111] orientations. The data were obtained using strain controlled 
loading in combination with an out-of-phase temperature cycle, for which the 
maximum compressive load occurs at the maximum temperature and maximum tensile 
load at the minimum temperature (Cycle I TMF). Although both the applied 
strain load and temperature had a balanced cycle with respect to time, an 
unbalanced hysteresis loop resulted because the material properties change 
with temperature. Since the Young's Modulus changes with temperature, that 
portion of total strain range which is inelastic was difficult to determine 
accurately. Although the final correlation for the tests used inelastic strain 
estimated from isothermal data at cycle mean temperature (Reference 14), the 
applicability of this to more general and unbalanced cycles needs to be 
evaluated. The inelastic strains could also be determined by analysis, 
provided accurate constitutive models can be developed. 

The Stress-Strain-Time (SST) model is identical in content to the Hysteretic 
Energy model except the loading terms are normalized to material strength 
capabilities (Reference 15). 


“f = A ' (%) b ' (■?)' ' T " (c - 7 > 


where Ae 

= 

strain range, 

e f 

= 

true strain at fracture. 

CT T 

= 

peak tensile stress, 

a y 

= 

yield stress. 



/Tq + Tj + T c f or T t > Tq 

T 


( T 0 for r T < r c 

T o 

= 

time per cycle spent in cycling, 

t t 

= 

time per cycle spent in tension hold, and 

T c 

= 

time per cycle spent in compression hold. 



The SST model was intended to provide a general framework into which the 
behavior of many different materials and load conditions could be 
incorporated. This is accomplished by allowing a choice of inelastic, elastic 
or total strain ranges for the correlating load parameter, depending on which 
strain provides the best correlation for the material and condition evaluated. 
Additionally, the stress term is used only when mean stresses are present, and 
frequency separation methods are permitted, although not recommended, under 
some conditions. 

C.3 Cumulative Damage Models 

The cumulative damage approaches assume that the plastic and creep components 
of inelastic strain cause damage which can be explicitly predicted and which 
define the state of the material. Damage is considered to be zero in the 
initial undamaged state and failure occurs when a critical level or limit is 
reached due to plastic, cyclic creep or monotonic creep deformations. A number 
of different definitions of damage and approaches for counting damage are 
used. Nearly all of these require use of a constitutive model to determine the 
portion of damage caused by a particular load condition. Life prediction 
models which use this approach are: Linear Time and Cycle Fraction, Ductility 
Exhaustion, Continuous Damage, Strain Range Partitioning, and Cyclic Damage 
Accumulation. 

C.3.1 Linear Time and Cycle Fraction Method 

One of the most widely used and highly developed cumulative damage life 
prediction schemes is the linear time, cycle fraction damage summation 
(Reference 16). A number of versions of this method have been employed. One of 
the most popular is that of code case 1592 of the ASME Boiler and Pressure 
Vessel Code (Reference 17). The basic equations for this method are shown 
below. 


where 


n 



= number of applied cycles at loading condition. 


Nd = life in cycles at loading condition, 

t = hold time at load and temperature, 


td = stress-rupture life at load and temperature, and 
D = allowable damage before failure. 


(C-8) 
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In some versions of the linear time and cycle fraction method, it was assumed 
that failure occurs when damage equals 1.0. In the ASME Pressure Vessel and 
Boiler Code, this assumption is dropped and the damage at failure is a 
nonlinear function of creep/fatigue life fractions. 

In order to use this method, correlations for predicting both Nj and t^ 
must be available for load ranges of interest. These usually are stress-time 
to rupture and strain-life relations similar to the Coffin-Manson model 
described above. Both are usually temperature dependent. 

For thermomechanical fatigue (TMF), the temperature is not constant and the 
prediction of and t^ is complicated. An effective method is not yet 
available for predicting for thermomechanical conditions, and, therefore, 
data for the temperature which yields a minimum life at a given strain usually 
are chosen when TMF data are unavailable. For the creep component of damage, 
however, the variable temperature can be dealt with using an integral equation 
(References 16 and 18) shown below. 

o 

where t r ( |S| ,T) = stress-rupture life at absolute stress level and 

temperature. 

Using these two methods of dealing with variable temperature, several 
different thermomechanical fatigue results have been successfully predicted 
(References 19 and 20). Isothermal results involving hold times also have been 
predicted (References 21 and 22). 

In spite of these successes, there is considerable evidence to suggest that 
the linear time and cycle fraction method is not satisfactory (References 10, 
22 and 23). Manson discussed some reasons for this at a 1978 AGARD conference 
(Reference 23). These pertain to describing the effects of complex damage 
histories and assessing damage due to compressive creep. The basic assumption 
for the linear model is that cyclic fatigue and creep damage components are 
mutually independent and can be calculated separately (Reference 24). This 
appears a consequence to the early observation that fatigue failures usually 
result in transgranular cracking while stress rupture cracks occur in grain 
boundaries in multigrained materials, and generally is not an accurate 
assessment of high temperature damage processes. 

Finally, since the linear model does not account for the different types of 
unbalanced creep, modeling effects of the observed tension/compression 
asymmetry in single crystals would clearly make the use of this model even 
more difficult. In principle, the directionality of fatigue and tensile 
properties could be accounted for in this method. But use of this model when 
the principal fatigue stress directions are not coincident with the creep 
directions is beyond demonstrated experience. 
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C.3.2 Ductility Exhaustion 

Ductility Exhaustion theory assumes that ductility is a material property and 
that fatigue and creep use up ductility to produce failure. In a given cycle, 
the ductility is debited a certain amount for the fatigue damage and a certain 
amount for the creep damage. The individual damage components are calculated 
using conventional fatigue life and creep curves and incrementally cumulated 
cycle by cycle. Failure results when the ductility has been exhausted. One 
version of this method used by Pratt & Whitney (Reference 24) is described as 
f ol 1 ows . 

The fatigue life is defined as a logarithmic relationship between total strain 
range and life: 


Acj • (N f ) 1/a = K (C-10) 

where Aej = total strain range, 



ef = 1/4 cycle intercept or fracture ductility, and 
a = material constant. 

The decrease in available ductility produced after n fatigue cycles is then 
written as: 

1 

a (C-ll) 


A decrease in available ductility due to creep is equal to the amount of creep. 

Ductility decrements from such cycling and creep can then be accumulated by 
stepwise progression along a creep- time curve to the point where the initially 
available ductility is exhausted. The creep curve is selected to be 
representative of the material behavior. The form of the equation for several 
current materials can be approximated by a power law. 

s = A t b (C-12) 



where t 


time, and 


A, B = material constants at stress and temperature. 
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In this method it is assumed the fatigue damage and creep damage are 
contributors to failure. As originally developed and currently implemented at 
Pratt & Whitney, thermal mechanical fatigue data are utilized to determine the 
fatigue ductility equation. Since a "faithful" test cycle could be used for 
obtaining these data, the cyclic creep strains are directly included in the 
fatigue data and the ductility exhaustion calculation accounts only for the 
effects of monotonic creep during hold. This method has been found quite 
satisfactory for certain cycles, especially those associated with turbine 
blades. The creep and fatigue damage interaction is' nonlinear and depends on 
the fracture ductility and the creep curve shape (Reference 24). 

In a recent investigation (Reference 3), the fatigue data were generated using 
rapid strain cycling resulting in practically no cyclic creep damage. For 
subsequent slower cycles and cycles including hold times, cyclic creep effects 
were added by the ductility exhaustion calculations and resulted in 
under-predicting life for these test conditions. Life trends observed in 
additional tests at varying R-ratios were found to be not predicted. These 
difficulties to account for the cycle shapes are not completely unexpected 
since: 1) the model has the same limitation to predicting compressive creep 
effects as the linear fraction method, and 2) no phenomenological correction 
is provided for maximum and mean stress effects. 

C.3.3 Continuous Damage 

The Continuous Damage prediction model (References 27 through 30) is a stress 
based model which assumes that the low cycle fatigue process can be described 
by a latent internal damage parameter (D ) which defines the progressive 
deterioration process taking place between the initial undamaged state and 
final failure. The deterioration model is that of the growth of assumed 
microdefects which are distributed throughout the material volume so that 
material strength can still be considered to be homogeneous on a macroscopic 
level. These microdefects may actually include effects of a number of 
metallurgical phenomena, including slip, ascent of dislocations, and inter- or 
intragranular microcracks (Reference 28). The failure is considered to be the 
initiation of larger cracks which affect the local stress distributions. 

A characteristic feature of the Continuous Damage model is the definition of 
damage growth in differential form for both the creep and fatigue damage 
components. 

dD = f ( 0 , D, T) dt + g (a f1 , F, D, T) dN (C-13) 

where D = damage parameter, 

f = creep damage rate, 

g = cyclic damage rate, 
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a - instantaneous stress, 

CTj,] = maximum stress for load cycle, 

a = average stress for load cycle, 

T = temperature, 

t = time for creep, and 

N = fatigue cycles. 

When the two damage rate terms, selected to define failure behavior of a 
particular material, and the loads are simply formulated, this expression can 
be integrated between suitable bounds and an analytical solution obtained. 
However, numerical solution techniques will be required when simple 
formulations are not possible. A nonlinear creep and fatigue life relationship 
results due to the presence of D in both parts of the differential equation. 

In the published work (References 27 through 30), the creep damage is modeled 
for constant or slowly varying loads using an equation based on those of 
References 25 and 26, modified to provide a more general nonlinear form 
dictated by the materials studied. Tensile and compressive anisotropy, if any, 
could be treated by this form (Reference 30). 


d0 c 

3T 


■[*]'■ 


(1-D) 


k(a) 


(C-14) 


where D c = creep damage, 

A+_ = material constants for tensile or compressive creep, and 

r, K(o) = material constants. 

Similarly, cyclic damage is defined (Reference 28). 


where a 
/3 
i-1 



material constants. 


(C-15) 
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These equations can be directly integrated between suitable bounds in order 
obtain life relationships for stress load conditions (References 25 through 
27). Applications to strain load conditions require use of constitutive models 
relating stresses and strains. Solutions for several special case load 
conditions are presented in Reference 30. The life relationships for strain 
controlled plastic low cycle fatigue and a saw tooth creep strain loads are 
identical to those previously described for phenomenological models in 
equations (C-l) and (C-2). 

For the general case of interacting creep and fatigue damage, differential 
equations (C-14) and ( C- 15 ) are substituted into equation (C-13) and solved by 
numerical methods. When creep and fatigue act simultaneously during each cycle 
(e.g., cyclic strain loads at high temperature or cyclic loads with strain 
controlled hold periods) the differential equation which must be solved can be 
rearranged as follows (Reference 31). 


dD 


(1-D) 


-k(a) 


[ 1 .( 1 - 


dN 


Ck ( a ) + 1] N c + ( 0 + 1) [1 - a ( 


a)] ‘ (1 - D) P N 


(C-16) 


where N c = cycles to failure under pure creep process, and 

Nf = cycles to failure due to pure fatigue. 

More recent work on Continuous Damage (References 31 and 32) discusses a 
general differential damage equation to replace equation (C-13). Use is made 
of two levels of damage, D a and Di to describe the condition of the 
material. D a allows consideration of fatigue damage "incubation" periods 
during which creep and fatigue damages are accrued separately, followed by a 
period of interaction. 


dD a = 9 ^ a M’ a ’ D » 


dD = f (a, D, T) • dt + H (1 - D g ) * g (a^, a, D, T) * dt (C-17) 

where H (1-D a ) =0 if D a < 1, 

= 1 if D a >_ 1. 

The Continuous Damage method was successfully applied to IN 100, several 
stainless steels, and OHFC copper. Using material data determined by pure 
creep and pure fatigue tests, highly nonlinear interaction effects could be 
predicted for symmetrical load and strain controlled conditions. An 
encouraging result was the ability to predict life trends when loading 
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sequence was changed. However, low frequency test results were consistently 
under-predicted in subsequent applications of this method to B1900+Hf and 
Hastelloy X (References 3 and 33), indicating further development may be 
necessary. Since the determination of material constants for this model is 
sensitive to the assumed constitutive behavior, constitutive modeling is one 
area where improvement appears necessary. 

Transient loads and temperatures, such as those which may occur during turbine 
engine acceleration and deceleration, have received only nominal consideration 
with very limited success in published literature. For such load conditions, 
damage definition is complicated because the differential damage expression is 
expected to require an extended form as follows: 

dO = U . dt + 7^ . da + ^ . dT ( C - 18 ) 

ot o a 61 

Equation (C-13), applicable to constant loads, includes only the time rate 
UD/at) portion of creep damage and a corresponding term for cyclic damage. 
Reference 35 recognizes this and proposes a modification to the differential 
creep damages in equations (C-13) and (C-14) by including a stress rate term. 
This approach nas not been sufficiently evaluated, however. 

dse of continuous damage analysis for modeling thermomechanical fatigue (TMF) 
life under slow cycles, i.e., negligible temperature rate ( 6D/<sT ) effects, has 
been attempted by some investigators (References 33 and 34). To account for 
varying temperature effects in model equations, an effective temperature is 
determined by considering the TMF cycle on an ultimate strength based 
nondimensional stress basis. The effective temperature is then used for life 
prediction. It was found that this method did not correctly predict Hastelloy 
X test results, although the same method had been successfully applied to IN 
100. The prediction difference was attributed to the relative strain rate 
sensitivities of the materials at their respective test conditions. 

C.3.4 Strain Range Partitioning (SRP) 

Strain Range Partitioning (References 23, and 36 through 39) is an extensively 
developed cumulative damage model which uses a strain range weighted linear 
damage cumulation method and results in a nonlinear life relationship that 
recognizes the differences between plastic and creep deformations, and between 
tension and compression loading. It had been previously noted that creep 
reduces cyclic life at elevated temperatures; however, considerable amounts of 
apparently conflicting evidence concerning the amount of life reduction were 
noted. This resulted in the conclusion that life reducing effects of creep 
largely depend on where in the cycle creep occurs, consistent with 
metal lograpnic findings that tensile and compressive cycle creep causes 
different types of failure, and evolved in developing SRP (Reference 36). 
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It is assumed that inelastic strain effects on life can be described by 
Coffin-Manson type relationships of the form described by equation (C-l), 
provided the hysteresis loop conforms to a specific form. Four characteristic 
forms are recognized and assumed to represent four basic components of damage. 
Life relationships for each are different and can be described using 
subscripts "p" and "c" to identify terms applicable to plastic and creep 
strain components respectively, as follows. 
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where 


pp = loop for high frequency cycling (creep = 0), 

cp = loop for tensile hold, 

pc = loop for compressive hold, 

cc = loop for compressive and tensile holds, and 

C and M = material constants. 


Loading cycles of most high temperature structures do not exactly conform to 
one of these four characteristic cycles. Fatigue life for these more complex 
cycles is calculated by partitioning the inelastic strain range into its four 
basic components. The partitioning method is defined for isothermal closed 
hysteresis loops for which the tensile and compressive inelastic strains are 
equal. In general, actual loading cycles will consist of three components, so 
that the sum of the components equals the inelastic strain range. Damage for 
each component is then added according to the interaction damage rule. 
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Strain range partitioning ( SRP ) has been successful in correlating a large 
amount of isothermal fatigue data involving unsymmetrical cycles and hold 
times for different materials, ranging from steels to various generations of 
nickel alloy (References 29, and 36 through 39). Although material and 
temperature independence for life relationships is not a requirement for the 
use of SRP, a method for collapsing data for various materials and 
temperatures is presented in Reference (40). This procedure provides the 
obvious advantages of minimizing total requirements, for fatigue data. 

Further refinements to SRP have been presented in two relatively recent 
documents (References 41 and 42). The SRP method as presented does not 
explicitly account for mean stress effects on life, and several investigators 
found difficulty in applying SRP to strong materials in the presence of mean 
stress (References 43 and 44). In some cases, it was necessary to use negative 
damage which resulted in mathematical difficulty. Reference 41 provides a 
method of explicitly correcting for effects of such mean stresses which can 
develop with unsymmetrical load cycles or material characteristics. 

A total strain range version of SRP is presented in Reference 42 in order to 
extend predictive capability to very low inelastic strains, prevalent in many 
practical fatigue problems where long lives are a requirement (e.g., turbine 
blades). Effectiveness of this version is strongly dependent on constitutive 
modeling accuracy. 

Use of SRP for thermomechanical fatigue (TMF) cycling may present 
difficulties, first, because of the temperature dependence of fatigue life for 
some materials (References 29 and 46) and, second, because the strain range 
partitioning process is not precisely defined for all TMF cycle hysteresis 
loops. The partitioning process requires that the tensile and compressive 
going inelastic strains are equal (Reference 36). This requirement is violated 
for hysteresis loops which are not closed because of mono tonic creep strains, 
such as those in a thermally cycled rotating turbine blade. Another case is 
the constant strain amplitude TMF cycle. 

Applications of SRP to TMF cycling where the basic strain/life relations are 
temperature independent have been presented for several materials in 
References 33, 39, and 46 through 48. In some cases (Reference 33), the 
unequaled compressive going inelastic (plastic) strain was assumed part of 
Ae pp ; in some others, the partitioning details are not explained. Generally, 
results of these applications were less than satisfactory. 

C.3.5 Cycle Damage Accumulation (CDA) 

Cycle damage accumulation is the newest cyclic life prediction method and 
currently still under development for uncoated B1900+Hf nickel base alloy 
(Reference 3). Cycle damage accumulation presents a new prediction approach 
which evolved as a result of evaluating various previous life prediction 
models for use on a typical high strength isotropic material for engine hot 
section components and is intended to provide improved creep-fatigue 
predictability. 
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Fundamentally, CDA follows concepts of both the ductility exhaustion 
(Reference 24) and continuous damage (References 27 through 30) models. 
However, the available ductility is redefined as primary creep in a material 
instead of the quarter cycle intercept. This evolved from metal lographic 
observations that dislocations which ultimately proceed to failure are created 
early in the load cycle and that the crack initiation process remains 
primarily transgranular, even though final rupture failure may be 
intergranular. The primary creep strain is considered indicative of grain 
deformation during initial creep stage. 

The solution technique integrates a damage rate over the loading history 
similar to that of the continuous damage method; however, a simplified damage 
rate definition is used. The approach is to determine damage at actual load 
conditions by applying damage ratios to previously tested reference 
conditions. This results in the following integral equation. 



where Tp = primary creep strain (ductility), 

dp = reference cyclic damage rate, 
dN R 

= A Ae-j- 


aj = maximum tensile stress, 

Ao = stress range, 

t = time, 

N = number of cycles, 

A, B, B', C' = material constants, and 
subscript R = reference condition. 



0 


( C-21) 


The initial development of this model was conducted using isothermal 
conditions at various frequencies, hold times and mean stresses. This resulted 
in a linear damage relationship as the damage rate is assumed constant through 
the cycle loading and the damage ratios are linearly added. It is not yet 
determined whether this relationship is a characteristic of the material 
tested, unlike most other materials which required nonlinear damage addition 
(References 24, and 27 through 30), or simply the status of currently 
incomplete model development. Similarly not determined is how the base metal 
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primary creep can be used to predict cracks in coated surfaces where the crack 
starts in the coating. Additional investigations to resolve these questions 
are currently in progress using variable temperature and TMF cycles and coated 
specimens. 

Although many model details and the general applicability of this model 
approach are still being determined, the use of the damage ratio concept 
provides the advantage of not having to evaluate absolute damage increments 
per cycle. In cycles that include significant amounts of time dependent 
damage, this is considered important due to the uncertainty associated with 
relating monotonic creep to the cyclic fatigue process (Reference 23). 

C.4 Crack Growth Models 


C.4.1 Damage Rate Model 

The Damage Rate Model (DRM) assumes that undamaged materials contain two types 
of microdefects, microcracks and cavities, and that fatigue can be modeled by 
crack propagation and cavity growth (Reference 49). The model is formulated so 
that both reversible and irreversible damage can be predicted: microcracks 
propagate under both tensile and compressive stresses but at different rates, 
while cavities will enlarge under tensile loads and shrink under compressive 
loads. The originally defined model was concerned with real transgranular and 
intergranular cracks, and grain boundary cavities, but actual applications 
could work with corresponding effective defects, similar to the continuous 
damage approach (References 27 through 30). This formulation can account for 
tensile and compressive damage differences. 

The growth laws are defined as follows: 
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for tensile stresses 

(C-22) 

for compressive stresses 
for tensile stresses 

(C-23) 

for compressive stresses 


where a = microcrack length, 
c = cavity length, 
e i pel = inelastic strain, 
einel = inelastic strain rate, and 
B, G, m, K, K c = material constants. 
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Life prediction equations result from integrating these rate equations over 
the specimen life from the initial microdefect sizes to the critical defect 
sizes when a macrocrack is formed, and using a microcrack and cavity growth 
interaction equation. Equation (C-24) defines cyclic life for an ideal cycle 
(Reference 11) under continuous cycling at constant frequency or strain rate. 
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c + T 

where A = ^ /In (af/a Q ) » 

a 0 = initial microcrack size, and 
af = critical microcrack size. 

The ideal cycle has been defined as having: 

peak tensile stress = peak compressive stress 

peak tensile inelastic strain = peak compressive inelastic strain, and 
e-j n ei = constant. 

When hold times are present, a linear damage cumulation method is used. 




( C-25) 


where Nf 1 = cycle life with hold times, and 

D , D , D , D = tensile and compressive creep damage terms for 
a T a c C T c c cracks and cavities (Reference 11). 

DRM was applied to two nickel base superalloys, Rene' 95 (Reference 11) and 
B1900+Hf (Reference 3), at isothermal conditions but did not result in 
satisfactory solutions. Applications to coatings or coated surfaces under TMF 
load conditions are not discussed in the surveyed literature. 


C.4.2 Fracture Mechanics 

Fracture mechanics analysis is generally considered to be an approach to 
describe life remaining in a part when a macroscopic crack or defect is 
present (Reference 50). The prediction method usually practiced is that of 
Linear Elastic Fracture Mechanics (LEFM) which involves describing local loads 
at the crack tip using stress or strain intensity determined from load 
distributions on the uncracked part and the crack size and geometry. Also 
required is the crack propagation law describing crack growth relative to the 
stress or strain intensity, as follows: 
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where aK = stress or strain intensity. 

For most materials, the propagation law includes a logarithmic linear relation 
for intermediate levels of aK, which can be readily integrated between the 
initial crack size and failure. However, for turbine blade materials and very 
long lives, the propagation law must be extended to low aK values, where the 
relation is nonlinear and numerical integration is required. Cook and Rau 
(Reference 51) reviewed the use of anisotropic fracture mechanics and arrived 
at a number of conclusions concerning the development of a life prediction 
method for anisotropic materials. They suggested that a reasonable life 
prediction model could be developed by utilizing the isotropic stress 
intensity factor(s) and empirically lumping all of the anisotropic effects 
into the material's response over the range of grain orientations. 

The success of using fracture mechanics for life predictions at elevated 
temperatures is directly related to knowing the crack growth laws for 
appropriate load conditions. Limited growth property tests have been reported 
for vane and blade materials for various temperatures and load cycles 
(References 52 through 54). 

Rau, Gemma and Leverant (Reference 52) found that crack propagation rates for 
nickel and cobalt base superalloys depend on mean stress, load cycle shape and 
temperature level. In another paper, Gemma, Langer and Leverant (Reference 53) 
investigated the effect of orientation on the crack propagation rate for the 
directionally-solidified nickel-base superalloy MAR M 200 subjected to thermal 
mechanical fatigue. They found that the crack propagation rate was strongly 
dependent upon the angle between the axis of loading and the grain growth 
direction. In addition, they showed that the orientation dependence on crack 
growth rate could be correlated using a normalized ratio of elastic moduli. 

The effects of stress dwells and varying mean strain on crack growth rates 
during thermal mechanical fatigue also have been investigated by Gemma, 

Ashland and Masci (Reference 54) for a conventionally cast and 
directionally-solidified nickel-base superalloy, MAR M 200. They showed that 
baseline primary creep rates correlated well with cumulative creep rates. They 
found no correlation of cumulative creep rates with secondary creep rates. In 
a NASA study, Bizon, Dreshfield and Calfo investigated the thermal fatigue 
resistance of a coated and uncoated directionally- solidified MAR M 247 
superalloy using a MACH 1 burner facility which permitted the leading o edge of 
the specimen to be thermally cycled between room temperature and 1070°C 
(1926 F) , Reference 55. They also found that the thermal fatigue resistance 
was strongly dependent on the grain orientation, with specimens having higher 
elastic modulus in the load direction showing correspondingly lower thermal 
fatigue resistance. Their data were in agreement with that obtained by Gemma 
et al (Reference 53) on MAR M 200 using crack growth rates under TMF 
conditions. 
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At least three applications of LEFM to coated turbine materials have been 
considered. The first of these involves directionally-sol idified turbine blade 
cracks which develop in the coating early during engine operation and the life 
is practically all propagation time (Reference 56). Measured crack depths were 
correlated with predictions. The second application considers actual defects 
in coatings on laboratory specimens (Reference 57). The defects which were 
characterized after specimen failure are treated as cracks which propagate 
initially in the coating and proceed into base metal. Local strain loads were 
calculated in the coating and used in the initial propagation phase. 

The third application also is for specimen analysis. Fifty percent load drop 
life is predicted in TMF specimens, assuming that coating cracks start early 
(Reference 58). Time for crack propagation in the base metal is then 
correlated with specimen life, using strain intensity calculation for observed 
crack shapes. A very significant result of this study was the correlation of 
lives for single crystal materials using an approach based on slip behavior of 
face centered cubic crystals. It was assumed that slip was equally distributed 
among the available slip directions. 

When significant macroscopic inelastic strains are present, LEFM is not 
applicable and the J-integral, C* or other nonlinear parameters must be used 
in place of aK (References 59 and 60). This has not been found necessary for 
turbine blade materials in which macroscopic inelastic strains do not occur 
under normal conditions, but is a requirement for weaker materials such as 
those used in combustor liners. Meyers (Reference 61) investigated several 
fracture mechanics analysis methods for predicting crack propagation life for 
Hasten oy X subjected to inelastic strains. The parameter which showed the 
most effectiveness in correlating high temperature and variable temperature 
crack growth data was Crack Opening Displacement calculated from the 
J-integral . 
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C.5 Metal lographic Investigations 

Many metal lographic studies have been published which shed light on the 
mechanisms of fatigue. Some studies concentrate on the 'larger scale' data, 
such as crack initiation locations, crack orientations, and modes of crack 
initiation and propagation; while other studies probe the 'finer scale' data 
of dislocation motion, the development of crystallographic slip bands and the 
subsequent nucleation of cracks. 

C.5.1 Crack Initiation Sites 

Many investigators have observed that crack initiation during fatigue is 
related to the movement of dislocations, and models have been proposed to 
explain this dislocation behavior and different types for slip in pure metals 
(Reference 62). Several investigators recently have studied relationships 
between slip and crack nucleation in model superalloys at room temperature 
(References 63 and 64). Crack nucleation and slip were observed at the surface 
of a smooth wall test specimen in one experiment (Reference 63). 
Crystallographic slip bands developed on the surface early during loading. 
Shallow microcracks were found in the slip bands even after a few cycles. 

After subsequent cycling, the slip line cracks became deeper and grew due to 
both extension and linking. A second type of microcrack also developed at an 
angle to the slip and normal to the applied stress load. This crack grew along 
a zigzag path due to extension only. 

Fatigue crack initiation in hot section anistropic superalloys (both single 
crystal and columnar grain) has been almost always associated with defects in 
the microstructure. These may be brittle phases, internal porosity, surface 
defects or environmental attack sites (References 13, and 65 through 72). The 
cracks initiate at the defects and propagate to failure under continued 
loading. Increased defect sizes cause shorter lives (References 65, 67 and 
68), but determining which particular defect initiates the crack depends on 
the material and on loading conditions. 

These results made use of uncoated specimen tests, similar to those used to 
obtain data for life prediction correlations. When tests were conducted on 
coated specimens, the cracking origin generally moved to the coating 
(References 13, 63, 66, 67 and 73). The limited published data on actual hot 
section cracking (References 13 and 56) and Pratt & Whitney experience also 
indicate that, in components using advanced anisotropic materials, cracks 
initiate in the coating. However, coated conventional multi grained (isotropic) 
materials may initiate cracks either in the coating or in the grain boundaries 
underneath the coating (Reference 68). 

The following table presents an overview of the initiation sites for the 
various materials which were included in this survey. 
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Crack Initiation Site Summary 


Material 


Initiation Site 

Carbides Porosity Oxidation Coati n 


brain 

Boundary Reference 



Low C 

Mar-M-200 DS 

PWA 1430 SC 

Coated PWA 
1480 SC 

Rene' 80 

Rene' 120, 

150 DS 

Rene' W4 SC 

Coated Rene' 
i\l4 SC 


X X 68-69 

X 13 

X 13, 73 

X 81 

XX 70 

X X 71-72, 77 

X 71, 77 


DS = Directionally Solidified 
SC = Single Crystal 
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C.5.2 Initiated Crack Orientations 


3oth crystallographic and noncrystall ographic transgranular cracks occur due 
to fatigue in anisotropic superalloys (References 6b and 68). The 
crystallographic cracks (sometimes termed stage I cracks) tend to form at low 
temperatures and high strain rates and are directly related to both material 
orientation and load orientation. The noncrystall ographic cracks (termed stage 
II cracks) are usually associated with high temperatures and low strain rates. 
Stage II crack direction is perpendicular to the applied stress direction but 
independent of material orientation. 

Planar bands of dislocations are generated at discontinuities (defects or 
previously initiated cracks) at all temperatures and strain rates (References 
65, 68 and 69). At the high temperatures/low strain rates, these dislocation 
bands disperse before stage I cracks can form, and subsequently stage II 
cracks are nucleated in a homogeneous dislocation array. The initiated stage 
II cracks propagate and subsequently transition to stage I cracking, depending 
on loading conditions (Reference 6y). However, at low temperatures/high strain 
rates, stage I cracks initiate within the slip bands before dislocations can 
disperse. 

In room temperature tests of single crystal Mar-M-200, crack initiation sites 
nearly always were MC type carbides oriented normal to the stress axis. These 
cracks formed when a tensile stress was applied. As the stress was increased, 
localized crystallographic slip occurred in the metal at both ends of the 
carbides and was followed by cracking after repeated loading. Similar slip and 
cracxing were also observed at a limited number of micropores present, but the 
cracking from micropores tended to be less severe (Reference 65). At higher 
temperatures, cracx initiation generally continued at both carbides and 
micropores; however, there was a greater tendency for cracking at the 
micropores (Reference 68). At the higher temperatures, slip became more 
homogeneous (References 68 and 74) and the initiated cracks were perpendicular 
to the load axis, particularly at slow strain rates. The cracks that led to 
failure continued to initiate internally at carbides and micropores in 
directionally solidified Mar-M-200 even when tested with an aluminide coating 
(References 65 and 68), although extensive cracking initiated in the coating 
and propagated a short distance into the base alloy. The coating cracks 
initiated and propagated normal to the stress axis again. 

In uncoated advanced single crystal superalloys, such as PWA 1480 and Rene' 

N4, fatigue cracks also initiate at defects (References 13, 71, 72, and 75). 
However, while the crack initiation in PWA 1480 occurred at subsurface or 
near-surface microporosity, the initiated cracks involved some amount of stage 
II type cracking adjacent to the microporosity at all temperatures (Reference 
13). The initiated cracks were perpendicular to the load direction independent 
of the material orientation for the <10Q> and <11 1> orientations tested. But 
the number of cycles to crack strongly depended on the material orientation. 
Similar results were reported for Rene' N4 (References 71 and 72) and an 
advanced model superalloy (Reference 76). At both 1400°F and 1800°F test 
temperatures, all six tested orientations of Rene' N4 involved a small region 
of stage II cracking associated with the crack origin, except the <11 1 > 
orientation at 140U°F, which showed an area of stage I planar cracking or cube 
slip planes. Also, Rene' N4 crack initiation at high temperatures occurred at 
the surface and was environmentally driven. 
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Cracks in both PWA 1480 and Rene 1 N4 initiated in the coating and were 
perpendicular to the applied load direction (References 13, 71, 72, and 77) 
when these alloys were tested in the coated condition. 

These results were previously observed in directionally solidified Mar-M-247 
material tested using wedge shaped test specimens in a burner rig (Reference 
55). Both coated and uncoated specimens cracked perpendicular to the specimen 
edge independent of the material orientation in the .specimen, which ranged 
from 0 to 30 degrees. However, the number of cycles to crack depended on the 
base alloy orientation; i.e., cycles to crack decreased as the angle between 
the specimen axis and the <001> direction increased. 

Environmental effects on fatigue cracking of superalloys are important in 
superalloys, for not only does surface oxidation present detrimental effects, 
but it can also provide advantageous effects on life (References 64, 67 and 
68). Experiments at moderate temperatures have shown that surface slip is more 
uniform and surface shear displacements smaller because of a tnin oxide film 
formation on surface slip offsets which leads to greater surface work 
hardening. At the higher temperatures, of course, the thin oxide benefits are 
overcome by local penetration of oxidation pitting which can act as crack 
initiation sites (References 68, 77 and 78). 

C.5.3 Coating Effects on Cracking 

Oxidation retardant coatings are used on hot section components mainly to 
provide protection against oxidation or corrosion during high temperature 
operation. However, the coatings also significantly influence the formation of 
cracks in the coated base alloy system and usually change the crack initiation 
site from the base alloy to the coating (References 68 and 79). 

The presence of the coating causes a bimetallic structure on the coated 
component surface which changes the surface deformation characteristics. The 
stresses developed in the coating can differ from those in the base alloy 
because of both constitutive behavior differences and thermal expansion 
property differences (References 56, 57, 79 and 80). 

Several investigators have considered the mechanical aspects of coatings. 
Reference 57 reported propagation of cracks in the coating and into the base 
alloy, resulting from thermal mechanical loads mainly due to coating and base 
alloy property differences. The cracks had initiated at coating defects 
believed inherent to the coating application method used. An experimental and 
analytical study (Reference 80) evaluated stresses in several coatings on an 
Armco iron bar under tensile load. The stresses caused coating failure along 
the coating grain boundaries. In another study (Reference 79), effects of 
strain resulting from coating application and subsequent service loads were 
considered. 

Other effects of a coating, previously mentioned, are eliminating environment 
related surface pitting from which cracks can initiate and crack initiation 
becoming controlled by coating instead of base metal fatigue properties 
(Reference 68). 
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Two distinct coating cracking modes have been observed, which indicate that 
two failure mechanisms may be important (References 13 and 71). Both modes 
result in cracks normal to the applied stress. The first mode is characterized 
by long cracks, completely or partially around the cylindrical test specimen 
circumference, whicn were found to initiate in aluminide coatings in 
isothermal fatigue at low temperatures and in-phase thermal mechanical fatigue 
at low temperature. These cracks are similar in appearance to limited 
ductility fractures and exhibited little evidence of inelastic deformation in 
the coating. At higher temperatures, a second mode is observed where both the 
aluminide and overlay coatings initiated small cracks locally. These cracks 
are thumonail -shaped and are preceded by coating deformation. 

C.5.4 Crack Initiation Criteria 

A wide range of cracking life criteria have been used in past investigations 
depending on material tested, test purpose and conditions. Specimen data, used 
in previous cracking life correlations, nearly always represented time to 
failure or impending failure indicated by a load drop. This circumstance was 
dictated at least in part by the internally initiated cracking which caused 
specimen failure but which could not be detected without destructive 
inspections. Experience shows, however, that a substantial amount of time is 
expended between crack initiation and final failure, controlled by crack 
propagation (References 13, 68, 71 and 75). Since including the propagation 
time with initiation life can lead to significant error, several investigators 
used periodic test interruption and replication to document surface conditions 
of specimens. When failure occurs from an externally initiated crack, 
initiation time for the critical crack can then be determined from the 
replicas (References 3, 13 and 77). 
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APPENDIX D 


TEST SEQUENCE OF CYCLIC CONSTITUTIVE TESTS 


The sequence of load for each test is listed in this appendix. The orientation 
of each specimen is also given in terms of deviation from the nominal crystal 
orientation and in terms of their Eulerian angles e and ip which are defined in 
Figure 3-7. 


Specimen Number = JA61 
Nominal Orientation = <100> 

Deviation from Nominal Orientation =6.7° 
Actual Orientation e = 6.50°, 4> = 0.68° 

Test Temperature = 427°C (800°F) 

Test Conditions 


Order of 

Strain Rate 

Nominal Strain Limits 

Number of Creep and 

Test 

( in/ in/ sec) 

( in/ in) 

Relaxation Test Conditions 

1 

io- 3 

+ .006 


2 

10" 3 

1 .008 


3 

10- 3 

1 .009 


4 

io- 3 

1 .010 
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Specimen Number = KA27 
Nominal Orientation = <110> 

Deviation from Nominal Orientation = 1.5 s 
Actual Orientation e = 43.60°, ^=0.50° 

Test Temperature = 427°C (800°F) 

Test Conditions 


Order of 

Strain Rate 

Nomi nal 

Strain Limits 

Number of Creep and 

Test 

( in/ in/ sec) 

( in/ in) 

Relaxation Test Conditions 

1 

10-3 

+ 

.003 


2 

lO- 3 

+ 

.004 


3 

10-3 

+ 

.0045 


4 

lO- 3 

+ 

.005 


5 

io - 4 

+ 

.005 


6 

io - 4 

+ 

.0053 
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Specimen Humber = LA66 



nominal Orientation = <111> 



Deviation 

from Nominal Orientation = 2.8' 


Actual Orientation e = 42.56°, 

ip= 33.28“ 


Test Temperature = 427°C (800°F) 



Test Conditions 



Order of 

Strain Rate Nominal Strain Limits 

Number of Creep and 

Test 

{ in/in/sec) 

(in/ in) 

Relaxation Test Conditions 

1 

10- 3 

+ .003 


2 

lO- 3 

* .0035 


3 

10-3 

+ .004 


4 

10-3 

+, .0045 


5 

lO- 3 

+ .005 

| 

6 

lO- 3 

+ .0053 


7 

10-3 

1 .0056 

| 

8 

10-3 

± .0061 


9 

lO" 3 

^ .0070 


10 

lO’ 3 

± .004 


11 

10-3 

+ .007 


12 

10-3 

1 -005 


13 

lO" 3 

i .007 


14 

lO" 3 

1 .006 


15 

10-3 

+ .005 


16 

10-3 

+ .0035 
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Specimen Number = MA27 
Nominal Orientation = <123> 

Deviation from Nominal Orientation = 1.2 s 
Actual Orientation e = 34.93°, ip = 15.40° 

Test Temperature = 427°C (800°F) 

Test Conditions 


Order of 
Test 

Strain Rate 
( in/in/sec) 

Nominal Strain Limits 
(in/in) 

Number of Creep and 
Relaxation Test Conditions 

1 

10- 3 

+ .003 


2 

10-3 

+ .004 


3 

10-3 

t .0045 

2 

4 

10-3 

1 .005 


5 

10-3 

+ .0053 

1 
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Specimen Number = JA44 
Nominal Orientation = <100> 

Deviation from Nominal Orientation = 5.7° 
Actual Orientation e = 4.89° , 4 >= 2.59° 

Test Temperature = 760°C (1400°F) 

Test Conditions 


Order of 

Strain Rate 

Nomi nal 

Strain Limits 

Number of Creep and 

Test 

( in/ in/ sec) 

( in/in) 

Relaxation Test Conditions 

1 

io- 3 

+ 

.006 


2 

IO" 3 

+ 

.007 


3 

io- 3 

+ 

.008 

5 

4 

io- 3 

+ 

.009 


5 

io- 3 

+ 

.010 


6 

io- 3 

+ 

.011 

5 

7 

5x10-5 

+ 

.011 


8 

io- 3 

+ 

.013 
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Specimen Number = KA26 

Nominal Orientation = <110> 

Deviation from Nominal Orientation =2.1° 
Actual Orientation e = 43.20°, ip = 1.03° 

Test Temperature = 760°C (1400°F) 

Test Conditions 


Order of 
Test 

Strain Rate 
(in/in/sec) 

1 

10- 3 

2 

lO- 3 

3 

10-3 

4 

10-3 

5 

10-3 

6 

lO" 4 

7 

lO" 3 

8 

lO" 3 

9 

lO" 3 


Nominal Strain Limits 
( in/ in) 

+ .003 

+ .003 

+ .004 

+ .005 

+ .006 

i .006 

± .006 

* .007 

+ .008 


Number of Creep and 
Relaxation Test Conditions 


10 
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Specimen Number = LA63 
Nominal Orientation = <111> 

Deviation from Nominal Orientation = 5.0° 
Actual Orientation e = 41.36", 4 >= 31.31" 

Test Temperature = 760°C (1400°F) 

Test Conditions 


Order of 

Strain Rate 

Nominal Strain Limits 

Number of Creep and 

Test 

( in/ in/ sec) 

(in/ in) 

Relaxation Test Conditions 

1 

10- 3 

+ .004 

10 

2 

io- 4 

+ .004 


3 

10-2 

+ .004 


4 

10-3 

1 .004 


5 

10-3 

+ .005 


6 

IO- 3 

1 .006 

7 
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Specimen Number = LA67 
Nominal Orientation = <111> 

Deviation from Nominal Orientation = 3.2° 
Actual Orientation e = 42.21“, 4 > = 33.00° 

Test Temperature = 760°C (1400°F) 

Test Conditions 


Order of 

Strain Rate 

Nominal Strain Limits 

Number of Creep and 

Test 

( in/ in/ sec) 

(in/in) 

Relaxation Test Conditions 

1 

lO- 3 

* .004 

Buckled on Loading 
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Specimen Number = MA25 
Nominal Orientation = <123> 

Deviation from Nominal Orientation = 2.0° 
Actual Orientation e = 35.75°, >P = 15.36° 

Test Temperature = 760°C (1400°F) 

Test Conditions 


Order of 
Test 

Strain Rate 
(in/ in/ sec) 

Nominal Strain Limits 
(in/ in) 

Number of Creep and 
Relaxation Test Conditions 

1 

10- 3 

+ .004 


2 

10-3 

LT> 

O 

o 

• 

+ i 


3 

10-3 

+ .006 

11 
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Specimen Number = JA66 
Nominal Orientation = <100> 

Deviation from Nominal Orientation * 5.3° 
Actual Orientation e = 3.90°, 4 1 = 2.79° 

Test Temperature = 982°C { 1800°F) 

Test Conditions 


rder of 

Strain Rate 

Nominal Strain Limits 

Number of Creep and 

Test 

(in/ in/ sec) 

(in/in) 

Relaxation Test Conditions 

1 

10-3 

+ .003 


2 

10- 3 

+ .004 

10 

3 

10-2 

* .004 


4 

lo-^ 

* .004 


5 

10- 5 

>ef 

O 

O 

. 

+ 1 


6 

lO- 3 

+ .004 


7 

10-3 

* .006 


8 

lO" 3 

+ .008 

12 

9 

lO" 2 

00 

o 

o 

• 

+ 1 


10 

*3* 

8 

o 

fl—fl 

i -008 

3 

11 

5xl0- 5 

CO 

o 

o 

• 

+ I 
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Specimen Number = KA22 
Nominal Orientation = <110> 

Deviation from Nominal Orientation = 5.4° 
Actual Orientation e = 41.55°, ip= 4.15° 

Test Temperature = 982°C (1300°F) 

Test Conditions 


Order of 
Test 

1 

2 

3 

4 

5 

6 
7 


Strain Rate 
( in/in/sec) 

Nominal Strain Limits 
( in/ in) 

Number of Creep and 
Relaxation Test Conditions 

10- 3 

+ 

.003 


10-3 

+ 

.004 

12 

10-2 

+ 

.004 


10-4 

+ 

.004 


10-3 

+ 

.005 


10-3 

+ 

.006 

10 

10-5 

+ 

.006 
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Specimen Number - LA68 
Nominal Orientation = <111> 

Deviation from Nominal Orientation = 1.4" 
Actual Orientation e = 43.97°, $ a 34.19° 

Test Temperature a 871°C (1600°F) 

Test Conditions 


Order of 
Test 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 


Strain Rate 
( in/ in/ sec) 

Nominal Strain Limits 
(in/ in) 

Number of Creep and 
Relaxation Test Conditions 

lO- 3 

+ .002 


lO- 3 

+ .003 

9 

10-3 

* 3 - 

o 

o 

0 

+ 1 


lO- 3 

+ .005 

9 

10-5 

+ .005 


10-4 

i .005 


10-3 

+ .005 


5x10-3 

in 

o 

o 

• 

+ ! 


lO" 3 

+ .005 


10-3 

* .006 

12 

lO" 3 

+ .010 
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Specimen Number = LA69 
Nominal Orientation = < 11 1> 

Deviation from Nominal Orientation = 2.0“ 
Actual Orientation 0 = 43.11°, ip= 34.02° 

Test Temperature = 982°C (1800°F) 

Test Conditions 


Order of 

Strain Rate 

Nominal Strain Limits 

Number of Creep and 

Test 

( in/ in/ sec) 

( in/ in) 

Relaxation Test Conditions 

1 

lO" 3 

+ .0015 

7 

2 

10" 3 

^ .002 

6 

3 

10- 3 

+ .003 

7 

4 

lO- 3 

+ .004 

12 

5 

10- 3 

+ .005 

9 

6 

10 - 3 

+ .004 


7 

lO" 3 

+ .°° 3 


8 

10" 3 

+ .002 

3 

9 

10' 3 

+ .0015 


10 

lO" 5 

+ .005 
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Specimen Number = LA64 

Nominal Orientation = < 11 1> 

Deviation from Nominal Orientation = 3.4° 
Actual Orientation e - 41.86°, <p = 33.07° 

Test Temperature = 982°C ( 1800°F) 

Test Conditions 


Order of 
Test 

Strain Rate 
(in/ in/ sec) 

Nominal Strain Limits 
(in/ in) 

Number of Creep and 
Relaxation Test Conditions 

1 

10“ 3 

+ .003 

7 

2 

10~ 4 

^ .003 

1 

3 

10“5 

+ .003 


4 

10- 4 

+ .003 

5 
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Specimen Number = MA23 
Nominal Orientation = <123> 

Deviation from Nominal Orientation = 5.1° 
Actual Orientation e = 35.60°, 4 > = 10.76° 

Test Temperature = 982°C (1800°F) 

Test Conditions 


rder of 

Strain Rate 

Nominal 

Strain Limits 

Number of Creep and 

Test 

(in/ in/ sec) 

( in/ in) 

Relaxation Test Conditions 

1 

io- 3 

+ 

.002 


2 

io- 3 

+ 

.003 

8 

3 

10-4 

+ 

.003 


4 

IO- 2 

+ 

.003 


5 

10-4 

+ 

.003 


6 

IO" 3 

+ 

.004 


7 

IO" 3 

+ 

.005 


8 

io- 5 

+ 

.005 


9 

IO- 3 

+ 

.005 

10 

10 

IO’ 2 

+ 

.005 


11 

IO' 3 

+ 

.0055 


12 

10-3 

+ 

.004 


13 

IO" 3 

+ 

.003 

3 

14 

5xl0“ 6 

+ 

.003 
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APPENDIX E 


SAMPLE DIGITAL FORMAT OF CYCLIC CONSTITUTIVE DATA 


Data from the cyclic constitutive tests are being recorded digitally as the 
testing proceeds so that virtually every cycle and hold period are preserved. 
The following is a sample of the format in which these data are being stored 
and will be reported. The data used in this example are taken from the 
stabilized response of specimen LA68, a <111> specimen tested at 871°C 
(1600°F). 

The format consists of two major "files", a HEADER file and a TRAILER file. 
The HEADER file gives a general description of the material and specimen 
geometry as well as the sequence of loading conditions imposed on the 
specimen. The TRAILER file gives the actual specimen response as 
stress-strain-time triplets for each "cycle" of each block of loading. The 
format has been designed so that keywords (e.g., MATL and SPECID) signal the 
type of data to follow. This format is more flexible and easier to read than 
strict card column formats. A more detailed description of the keywords 
follows. 


Keyword 


Description 


HEADER Signals the beginning of the HEADER file 

TITLE General description of test 

AN I SO 
THETA 
PSI 
GAMMA 

COMMENT Descriptive comment 

MATL Signals that general material data follows 

CODE Material specification 

HCODE Heat code 

HTSPEC Heat treatment specification 

GEOMETRY Signals that specimen geometry follows 

CONFIG Specimen part number 

GEOMUNIT Units used to describe the specimen dimensions 

OD Outer diameter 

GAGELN Gage length 

JOB Signals that project identification follows 

JOB ID Project identification 

SPECID Specimen identification (serial or specimen number) 

DATE Date of test 

BLOCi Signals that loading conditions for the ith block of loading 

will follow. A new BLOCi (sequentially numbered) is used each 
time the imposed loading conditions are changed. 


Signals that anisotropic material description follows 
Eulerian angles describing the specimen orientation 
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Keyword 


Description 


NCYC 


ISOTEMP 

TESTTYP 

UNITS 


STRN1UNIT 
STRS1UNIT 
TUNIT 
TIMUNIT 
CONVERT 
STRN1CON 
STRS1C0N 
TCON 
TIMCON 
TRAILER 
CYCLE 

BEGTIME 
STRESS 1 * * 
STRAIN 1 * * 
TIME * * 


The generic "cycle" number at which this loading condition 
begins. The "cycle" number corresponds to CYCLE in the TRAILER 
file describing the response under BLOC 1 loading conditions. 
Note that load or strain hold periods are also considered to 
be generic "cycles" in this format. 

Signals that the test was isothermal and is followed by the 
test temperature 

Signals that the test is a constitutive test 
Signals that the units for reporting the data follow. These 
units remain in effect for all tests and are specified only in 
the BLOC i description. 

Strain units 
Stress units 
Temperature units 
Time units 

These data are included to facilitate subsequent conversion of 
the data to other units. These conversion factors are specified 
only in the BL0C1 description. 

Signals the beginning of the TRAILER file 

The generic cycle number for the stress-strain-time triplets 

that follow 

The actual test time at the beginning of the cycle 
Signals that stress values will follow in free field format 
Signals that strain values will follow in free field format 
Signals that time values will follow in free field format 
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ORIGINAL PAS E IS 
PF POOR QUALITY 


HEADER 

..title cwsTirynvE_iEST _pwa 1*80 <iijl> nominal 

ANISO 


THETA 

= *3.97 

PSI 

a 39.19 

GAMMA 

a 0 

COMMENT 2 FROM < 111 > UNCOATED 

MATL 


CODE 

* PVU1A80 

HQODE 


HTSPEC 

= 1975F/9HR+1600F/32HR 

GEOMETRY 


CONFIG 

= LED *178* 

GEOMUNIT 

= INCHES 

00 

= .3 

GAGELN 

= 1 

JOS 


JOBID 

a NAS3-23929 

SPECID 

- U-66, T-729-8* 

OATE 

a 0B/30/8* 

BL0C1 


NCTC 

a 1 

ISOTEMP 

3 1600 

TESTYP 

a C 

COMMENT MINSTRAIN a -5.00000E-03 MAXSTRAINa 5.00000E-03 

COMMENT. STRAIN 

RATE_*_I . fiOOOOE-05 

UN JL 1 3 

STRN1UNIT 

a IN/IN 

STRS1UNIT 

a PSI 

TUNIT 

a F 

TINUNIT 

a SEC 

CONVERT 


STRN1C0N 

a 1.0 

STRS1C0N 

a 1.0 

TCON 

a 1.0 

TIMCON 

a 1.0 

BLOC 2 


NCTC 

a 2 

ISOTEMP 

a 1600 

TESTYP 

a C 

COMMENT MINSTRAIN * -5.00000E-03 MAXSTRAINa 5.00000E-03 

COMMENT STRAIN 

1 RATE a 1.00000E-0* 

BLOC 3 


NCTC 

a 3 

ISOTEMP 

= 1600 

TESTYP 

a C 

’ c OnM EN f n 1 NSTk A in = - 5 . 00 0 Ou E - 0 3 HAXSTRAIN= 5.00000E-03 

COMMENT STRAIN RATE = I.00000E-03 

BL0C4 

* 

NCYC 

a * 

IsBTeHP * 

a lSUO 

TESTYP 

a C 


COMMENT STRAIN HOLO AT 3.90B16E-03 
COMMENT 

BLOCS 
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NCYC s 5 

. _ ISOTEMP = 1600 __ 

TESTYP ~' = C ' ~ 7" — — 

COMMENT LOAD HOLD AT 15722.9 
COMMENT 

_CYCIE = 1 BEGTIME = 2097.9 

TRAILER ‘ 

STRESS 1 * * 

-96131.9 -98032.5 -99919.7 -51988.1 -52529.7 

-59252.5 -559§2_-57189.7_:5857.2 -59935.2 

-60990.9 -61682 -62718.7 -63755.3 -69792 
-65655.9 -66519.8 -67210.9 -67556.5 -68593.1 
-68938.7 -69802.6 -69975.9 -70839.3 -70839.3 

-71357\6 ;60990.9 -53093 ..1_-55L<>S -38356.9 

-31100.2 -29361.8 -17950.6 - ii 0885 .7" 9665. 02 
1382.23 6738.37 12099.5 16932.3 21597.3 

25571.2 29717.9 33000.7 36283.5 39566.3 

91985.2 99922.5 97 168 - 6 99587.5 51988.1 __ __ 

53561,9" 55962 "56899.2 58572 60127 

61509.2 62891.9 69100.9 69969.8 65828.7 

67038.1 67729.2 68765.9 69957 70320.9 

71012 71530.9 72098.7 72399.3 72912.6 

'" 73258.2 '79 i 2 2.1 "79299T8~ 79986 75 158/7 

75509.3 65828.7 58226.9 51192.5 99231.3 

37993 31100.2 29016.2 17623.9 11576.2 

5356.19 -395.557 -5879.98 -10712.3 -15895.6 __ __ 

- 1 986 97 6 -29016.2 -27299 -3023673 -3317375' 

-35592.9 -38189.1 -90089.7 -92330.8 -99058.6 
STRAIN 1 * » 

9.39668E-05 -1.66097E-09 -3.S1739E-09 -5.76953E-09 -7.66976E-09 __ 

-97 77039E-09 -1 .17733E-03 -T. 38251 E -03 -l'7'5'9766E-03 -1.79775E-03 
-2 . 00293E-03 -2.18857E-03 -2.90890E-03 -2.60870E-03 -2.81876E-03 
-3.01917E-03 -3.22923E-03 -3.92991E-03 -3.63959E-03 -3.89959E-03 
-9.09983E-03 -9.25989E-03 j;9.95530E-03 -9.66098E-03 ^9.86S66E^03 
-5. 051292 -03 -9.875'93E-03 -976751 3E- 0 3 '-979650 7E -03 -5.2S989E-03'"" 

-9 , 05971E-03 -3.85992E-03 -3.65913E-03 -3.95872E-03 -3.25359E-03 
-3.03859E-03 -2.83830E-03 -2.63312E-03 -2.93283E-03 -2.22276E-03 
-2.01759E-03 -1.807S2E-03 -1.61700E-03 -1.91182E-03 -1.19687E-03 
- i 7 00 1 971-03 -7 . 96 2 87 E- 09 -6'.05765E-09 -3.'9b'8i6E-09 -1 . 90523E-09 
2 .992602-05 2.19839E-09 9.20127E-09 6.25305E-09 8.90259E-09 

1 . 09055E-03 1 . 29573E-03 1.95579E-03 1.69631E-03 1.89172E-03 

2 . 06 155E-03 2.26185E-03 2.98168E-03 2.67709E-03 2.88227E-03 

3.'0?2^E-03' 37282851-03" 3”. 5(55691-1)3 "3.5'$32i'E-53 '375oSl7Po3 

9. 10357E-03 9 . 31 363E-03 9.51392E-03 9.71910E-03 9.91939E-03 

5. 12957E-03 9.99382E-03 9.79353E-03 9.53835E-03 9.33317E-03 

9. 12799E-03 3.92281E-03 3.71275E-03 3.51296E-03 3.32193E-03 

_ r. 1 62 10 e - 6 3 — 2'7 592 651 -15 3 — 5 - 15 J’— "2 V35 12 2E -0'3i — 2’ ;i27I'S'2 

2 . 09086E-03 1.88569E-03 1.67562E-03 1.98510E-03 1.27509E-03 

1 . 069S6E-03 8.69680E-09 6.59502E-09 9.59323E-09 2.58915E-09 

TIME * « 

““ T 25 .979 90.9^5 'ICXzr 81.506“ 

102.37 122.899 193.318 163.792 189.266 

209.79 225.219 295.688 266.162 286.636 

307.11 327.589 398.058 368.532 389.006 

— 509753 529 '.'955 55 07528 5707902 59T7376 
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OF POOR QUALITY 


511.85 

532.324 

552.798 

573.272 

593.746 


614.22 

634.694 

655.168 

675.642 

696.116 


716.59 

737.064 

757.538 

778.012 

798.466 


818.96 

839.434 

859.908 

880.382 

900.856 


921.33 

941.804 

962.278 

982.752 

1003.23 


1023.7 

1044.17 

1064.65 

1085.12 

1105.6 


1126.07 

1146.54 

1167.02 

1187.49 

1207.97 


1228.44 

1248.91 

1269.39 

1289.86 

1310.34 


1330.81 

1351.28 

1371.76 

1392.23 

1412.71 


1433.16 

1453.65 

1474.13 

1494 . 6_ 

15.11L.0S_ 


1535.55 

1556.02 

1576.5 

1596.97 

1617.45 


1637.92 

1658.39 

1678.87 

1699.34 

1719.82 


1740.29 

1760.76 

1781.24 

1801.71 

1822.19 


1842.66 

1863.13 

1883.61 

1904.08 

1924.56 


1945.03 

1965.5 

1985.98 

2006.45 

2026.93 


CTCLE = 

2 

8EGTIHE * 409 

.2 


STRESS 1 

* K 





-51488.1 

-53907 - 

55980.3 - 

58053.6 - 

59781.4 



-61682 -6340975 -6496475 -6669275 -66247.6 


-69802.6 

-71357.6 

-72567 -73949.3 -75331.5 

-76368.2 

-77577.6 

-78614.3 -79823.8 

-80514.9 

-81551.5 

-82415.4 

-83279.3 -84143.2 

-84834.3 

-85525.5 

-76195.4 

-67210.9 -59263.1 

-51833.6 

-44749.7 

-37493 -30581.8 -23497.9 -16759.5 

-10366.7 

-3801.13 

2418.9 8120.6 I 

13822.3 

19005.7 

24016.2 

28681.3 33173.5 

37320.2 

40775.8 

44231.3 

46823 50105.6 52697.5 

55289.2 

57362.5 

59608.6 61854.8 

63409.8 

65310.3 

67038.1 

68420.4 70320.9 

71530.4 

73085.4 

74122.1 

75504.3 76368.2 

77577.6 

78959.9 

79996.5 

80687.6 81551.5 

82588.2 

83279.3 

83970.4 

84834.3 85525.5 

86043.8 

86562.1 

78441.5 

70148.1 62373.1 

55289.2 

48205.3 

41466.9 

34383 27471.8 

20906.2 

14340.6 

7947.82 

2073.34 -3973.91 

-9675.61 

-15204.5 

-20215.1 

-24880.1 -29199.6 

-33519.1 

-37147.4 

-40430.2 

-43540.2 -46477.5 

-49241.9 


STRAIN 1 * * __ 

7732786e- 05 - 1 . 270 15E-04 -3. 32 193E-04 -5.42257E-04 -7. 37665E-04 
-9.42843E-04 -1.14802E-03 -1.35320E-03 -1.55838E-03 -1.75379E-03 
-1.96873E-03 -2.17391E-03 -2.37909E-03 -2.57450E-03 -2.78456E-03 
-2 . 98974E-03 -3.20469E-03 -3.39033E-03 -3.60528E-03 -3.81045E-03 
^4 .61 563E -03 -4 . 2 1593T-03"* -574T599I -ti'S ^.S'S'lITE-fa -47525581^03“ 
-5.02198E-03 -4.91451E-03 -4.70933E-03 -4.49927E-03 -4.29897E-03 
-4.09868E-03 -3.89839E-03 -3.68344E-03 -3.48803E-03 -3.2828SE-03 
-3. 08256E-03 -2.86761E-03 -2.66732E-03 -2.46702E-03 -2.26673E-03 
■"-O-5'6T7?-5r^I7^5S5r^3^55ir5rr-7r"-17^'561I-0'3~r':2^5£ : 53 , ™ — 
-1.040S5E-03 -8.402S4E-04 -6.35076E-04 -4.29897E-04 -2.24719E-04 
-l . 95408E-05 1.75867E-04 3.85931E-04 6.00879E-04 8.01172E-04 

9.96S80E-04 1.21153E-03 1.42159E-03 1.60723E-03 1.80752E-03 

- T0I735E-5S "l .13255 E-5 3 27'4 2356E - o!~ I765S12E-CTI *'2.F3835E^3 
3.04348E-03 3.24866E-03 3.44895E-03 3.64436E-03 3.86419E-03 

4.06448E-03 4.2S989E-03 4.47484E-03 4.68002E-03 4.88520E-03 

5.08061E-03 4 . 98290E-03 4.78261E-03 4.57254E-03 4.37225E-03 

5"."I5?I'9r^3'’“T: i J'6'67ST-5T~1':75r6'0r:03' - T:5'45?5E-DT“ , 7732r5?6F= , 0‘3 
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3.15584E-03 2.95066E-03 2.74060E-03 2.54030E-03 2.33512E-03 
_i.»22477E ; 03_l jM 72fiA7E-0|„ i .52<kt8E-9,3 J. .21900E;^3 
1 . 10405E-03 9.08647E-04 7.08354E-04 57o3“i'75E-04 2 . 93lI"2E-04~ _ 
TIME » * 

0 2.046 4.092 6.138 6.184 


10.23 

12.276 

14.322 

16.368 

18.414 


20.46 

22 .506 

24.552 

26 .'598 

28.644 


30.69 

32.736 

34.782 

36.828 

38.674 


40.92 

42.966 

45.012 

47.058 

49.104 


JU5„ 

53.196 

55.242 

52-M 8_ 

-59*334 


61.38 

63.426 

65.472 

67.518 

69.564 


71.61 

73.656 

75.702 

77.746 

79.794 


81.84 

83.886 

85.932 

67.978 

90.024 


92.07 

94.116 

96.162 

98.208 

100.254 


102.3 

104.346 

106.392 108.438 110.484 



112. 

53 

114. 

576 

116 

.622 

118. 

,668 

120. 

,714 

122. 

,76 

124. 

,806 

126 

.652 

128. 

,898 

130. 

,944 

132. 

,99 135, 

,036 

137 

.082 

139. 

,126 

141. 

,174 

143. 

.22 

145, 

.266“ 

~147 

.312 

149. 

.358 151. 

,404' 

153. 

,45 

155, 

,496 

157 

.542 

159, 

.588 

161, 

.634 

163. 

,68 

165, 

.726 

167 

.772 

169 

.818 

171, 

.864 

173, 

.91 

175, 

.956 

178 

.002 

180, 

.048 

182, 

.094 

184, 

:'i4 

186 , 

.186 

188 

.232' 

190, 

.'278 

192, 

[324 

194, 

.37 

196 

.416 

198 

.462 

200 

.508 

202 

.554 

CYCLE * 

3 


BEGTIME = 

60 




STRESS_1 *_•_ _ 

-52351.9 -55980.3 -59608.6 -63064.2 -66174.2 
-68938.7 -71875.9 -74294.8 -76541 -78959.9 
-80860.4 -83106.5 -84834.3 -86389.3 -88289.9 
-89672.1 -91054.4 -92609.4 -93991.6 -95373.8 


-96410.5 -97792.7 -98829.4 -99693.3 -100730 
-101594 -102285 -95028.3 -84316 -75677.1 


-67556.5 

-59954.2 

-52697.5 

-45440.8 

-36184.1 

-31100.2 

-24361.8 

-17277.9 

-10539.5 -3801.13 

2591768 

8984749 

14686.2 20387.9 

26780.7 

31964.1 

36801.9 

41466.9 

45959.1 

50278.6 

54252.5 

57880.9 

61336.4 

64619.2 

67729.2 

70493.7 

73085.4 

75677.1 

77750.4 

79651 

d'1378.6 

83279. 3~ 

85007.1 

86562.1 

88117.1 

89326.6 

90708.8 

92091 93473.3 

94509.9 

95719.4 

96928.8 

98138.3 

98829.4 

99866.1 


100730 101767 96410.5 86216.6 77923.2 

-7049377 “5289174 "'SsTSoTS “4975971 ' 4Tf5s72? 

35246.9 28508.5 21597.3 15031.7 8466.16 

2073.34 -4146.69 -9848.39 -16414 -22115.7 
-27817.4 -32655.2 -37665.8 -42158 -46477.5 

■■fmxiirr*"* 

2.83341E-04 7.61632E-05 -1.41671E-04 -3.41964E-04 -5.47142E-04 

-7.47435E-04 -9.42843E-04 -1.15779E-03 -1.35809E-03 -.1 .55838E-03 
-1.76844E-03 -1.96873E-03 -2.17880E-03 -2.37909E-03 -2.57450E-03 
■-TT7 54551-71 -7.77974r : 73-^37l9773E^0T*-T:4fl7l7I-ir3 =7.77 s25e=73 
-3.80557E-03 -4.01075E-03 -4.2U04E-03 -4.40645E-03 -4.62140E-03 
-4.82656E-03 -S.02198E-03 -4.97802E-03 -4.73376E-03 -4.52369E-03 
-4.31851E-03 -4. 10845E-03 -3.91304E-03 -3.70787E-03 -3.50269E-03 
-=1TST75IE -C3 - 3T m nT-gT^gTgg715r=g3-=7T5g 6B6 E -03 
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WBWB Sa mm aaw 


* "10175E-Q3 2.60o7wt**°^ 

;• i :; 5 , {-. llml- » « :Uu*-.s 


9.89612E-03 

9.92S99E-03 


5.05618E-03 

9.22570E-03 


2 .388S6E-03 

2 . 17880E-03 

1.36786E-03 

1.16756E-03 

TIME * * 

n ? .6 i A — 

1 1.2 1.9 

1.6 1.8 

2 2.2 2.9 

2.6 2.8 

3 3.2 3.9 

3.6 3.6 

9 9.2 9.9 

9.6 9_..8 

“ T 5.2 5.9 

5.6 5.8 

6 6.2 6.9 

6.6 6.6 

7 7.2 7.9 

7.6 7.8 

8 8.2 8.9 

8.6 6.6 

" 9 9.2 9.9 

9.6 9.8 


5 06106E-Q3 9.89123E-03 9.69099E-03 

9. 02052E-03 3.81095E-03 3.60SUE-03 

? . 99963E-0.1 _£• 28945E i a3„2^.9^9E-0i_ 
1 .77821E-03 1.56615E-03 

7.52320E-09 5.97192E-09 


1.99805E-03 

9.67269E-09 


JO. 9 
11.9 


10.6 

11.6 


11.8 


12 12.2 12.9 12 . .6 . 12 . 8 — 

<— - t" i?* /. i i . a 



_BEGTIME-„; _283 i 2_ 


52179.2 52179.2 

52351.9 52351.9 

52006.9 52006^9. 

51833.6 51833.6 

51660.8 51660.8 

51986.1 51988.1 


52351.9 
52179.2 

52006.9 


52351.9 
52179.2 


52351.9 

52179.2 - 

5 183 3 .6 51833.6, 

- 5166<h8 51833.6 

51315.3 51968.1 

51988.1 51315.3 


"51660.8 
51315.3 

.«««««— 

50969* 7 50969.7 50796.9 50629.2 50796.9 

50969.7 so796<9 50529.2 

50629.2 50951.9 

'~£0278".6 50278 ■ 6 

50276.6 50278.6 

99933 50105.8 99933 

99933 50105. 8_ 99933 

"'50To'5T8 '9 9933“ 99933 

" UU . Tm . j "™*!?” 5 «»«.. 

-4 


50796.9 50629.2 

50951.9 50629.2 

5095173 50951“ 9 

50278.6 50278.6 

50278.6 50105.8 

99933 _99933_ 
99933 “99933 


50951.9 

50951.9 

“50951.9 

50105.8 
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49587.5 

49414.7 

49587.5 

49241.9 

49414.7 

4941.4. 7_ 

49414. 2_ 

49414. 7_ 

49414.7 

. 1.4.2... _ _• ...._ 

49414.7 

49414.7 

49414.7 

49414.7 

49241.9 

49069.2 

49241.9 

49414.7 

49241.9 

49414.7 

49069.2 

49069.2 

49069.2 

49069.2 

49069.2 

49069.2 

49069.2 

49069.2 

49241.9 

49069.2 

48896.4 

49069.2 

49069.2 

46896.4 

48896.4 

48396.4 

49069.2 

48896.4 

48896.4 

49069.2 

48396.4 

48896.4 

48896.4 

48896.4 

48896.4 

_48Z?3,.9. 

_49069,? 

-48896 , 

..48596 .,4 

.4.9.059., .2 

48896.4 

48896.4 

46896.4 

48896.4 

48896.4 

48723.6 

48723.6 

48896.4 

48723.6 

48723.6 

48896.4 

49069.2 

48723.6 

48723.6 

48723.6 

48723.6 

48550,5 

48723.6 48723.6 

48723.6 

48550.8 

48550.8 

48550.8 

48550.8 

48550.8 

48550.8 

48550.8 

48550.8 

48378 48550.8 

48550.8 

48550.8 

48550.8 

48550.8 

48550.8 


46550 .8 48723.6 48378 48378 48378 SJ .;_ _ _ 

48378 48550.8 48378 48550.8' '48378' ‘ 

48550.8 48378 48378 48378 48378 

48378 48378 48378 48550.8 48205.3 

48378 48378 48205. 3_48205. 3 48205^3 

48378 48378~ 48205 i 3 48205.3 48378 

48378 48205.3 48205.3 48378 48378 

48378 48205.3 48378 48205.3 48205.3 

_48205 . 3 48205^3 48205.3 48205.3 _48205.3 __ 

48205.3 '48378 48032.5" 4820573 ""48032.' 5' 

48205.3 48032.5 48032.5 48205.3 48032.5 

48205.3 48032.5 47859.7 48032.5 

STRAIN I » » _ _ 

3. 90816E-03 3.913091-03 3.91304E-03 3.91793E-03 3.91793E-03 

3.91304E-03 3.91793E-03 3.91304E-03 3.91304E-03 3.91793E-03 

3.92281E-03 3.92281E-03 3.91304E-03 3.92281E-03 3.91793E-03 

V.91304E-03 3.91304EM13 _3.91793E-03 3.90816E-03 J.91793E-03 

3.91304E-03 3.91793E-03 3.91793E-03 37913046-03 3.91793E-03"' 

3.91793E-03 3.9I304E-03 3.91793E-03 3.91304E-03 3.90816E-03 

3.91304E-03 3.91793E-03 3.91793E-03 3.92281E-03 3.91793E-03 

3.91304E-03 3. 91 793E-03_3. 91304E-03_3. 91304E-.03 3 1 91304E-03_ 

3791 3041-03 5T91 304E-03 3. 9 1 304 E -03 3'.'91304E-03 3. 9T793E-03 

3. 91793E-03 3.91793E-03 3.91304E-03 3.91304E-03 3.91793E-03 

3. 90816E-03 3.91304E-03 3.91793E-03 3.91304E-03 3.91304E-03 

3. 91 793E-03 3.91793E-03_ 3.91304E-03 _3.91793£-03_ 3.91304E-03 

3 . 9'1 304E-03' 3 . 9 1304E-03 3.9T304E-03 3 .'9225IE-03 '3.91304^03' 

3. 91 304E-03 3.90816E-03 3.90816E-03 3.90816E-03 3.91304E-03 

3. 91 793E-03 3.90816E-03 3.91793E-03 3.91304E-03 3.90327E-03 

3.91304E-03 3.90327E-03 3.91304E-03 3.91793E-03 3.90327E-03 

TT^'oa i 6 e'-'o 3 3"79T3'04E^'5'S 3T?oSi5S'-?5 3T5TSME-03 

3.90816E-03 3.91304E-03 3.90816E-03 3.91793E-03 3.91304E-03 

3.91 304E-03 3.90816E-03 3.90816E-03 3.91304E-03 3.90816E-03 

3.91304E-03 3.91304E-03 3.91304E-03 3.90816E-03 3.91304E-03 

' 3: 9i 304E-03 ”3 . 51793 E - 0 3~ 3 1"§ 1 SO 4 1 - 0 S " 379081SE-53 "'5. 373 2 "7 £^3 
3. 90316E-03 3.91793E-03 3.90327E-03 3.91304E-03 3.91304E-03 

3.91793E-03 3.91304E-03 3.92281E-03 3.91304E-03 3.91304E-03 

3.91793E-03 3.90816E-03 3.90327E-03 3.91304E-03 3.91793E-03 

— 3:''9Wrsr^3''^:'908I?E-0T"''3;'9r30'4E-0'3'" _ ir'9'n05JE-O3~7:WJ27E-0'3 
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3. 91793E-03 

3.90816E-03 

3.90816E-03 

3.90816E-03 

3.91304E-03 

3. 90816E-03 

3 . 90816E-03 

3.913042-03_ 

1*913046-0.3. 

_3-9.9816.ErS 3 . 

3 . 90816E-03 

3.91304E-03 

3. 91304E-03 

3. 91304E-03 

3 . 90816E-03 

3. 9Q816E-03 

3. 91304E-03 

3.90816E-03 

3.90816'E-03 

3.90816E-03 

3.90327E-03 

3.91304E-03 

3.91304E-03 

3 . 91304E-03 

3.91793E-03 

3 . 91793E-03 

3 . 91 304E-03 

3 . 91304E-03 

3. 90816E-03 

3 . 90816E-03 

3.91304E-03 

3.91304E-03 

3.91304E-03 

3.90616E-03 

3.91304E-03 

3 . 91 304E-03 

3 . 91304E-03 

3 . 91793E-03 

3 . 91793E-03 

3. 91304E-03 

3. 91 304E-03 

3. 91304E-03 

3.91793E-03 

3.91793E-03 

3. 91 793E-03 

.3-2i3p9E-03_ 

_3.,aiZ235.?.oi_ 

1*9X3.041- 9.3. _1 • ,9128.VEz93. 

...3*3X1946-03 

3.92281E-03 

3.92281E-03 

3. 91304E-03 

3.90816E-03 

3.91793E-03 

3.92281E-03 

3.91304E-03 

3. 91793E-03 

3.91304E-03 

3.91793E-03 

3.91304E-03 

3.91793E-03 

3.91793E-03 

3.92281E-03 

3.92281E-03 

3.91793E-03 

3,91793E-Q3_3,9 1 793E - 03... 

3.91793E-03 

3.92281E-03 

3.92281E-03 

3.92281E-03 

3 . 91793E-03 

3.91793E-03 

3. 90816E-03 

3.91793E-03 

3.92281E-03 

3.91793E-03 

3.90816E-03 

3.91793E-03 

3. 91 304E-03 

3. 92281E-03 

3. 91793E-03 

3.91793E-03 

3.92281E-03 

3.91304E-03 

3. 91 304E-03 

3.91793E-03 

3. 91793E-03 

3.91304E-03 

3.91793E-03 

S.91793E-03 

3.91304E-03 

3.9V304E-03 

3.92281E-03 

3. 91793E-03 

3.91304E-03 

3.91304E-03 

3. 92770E-03 

3. 92281E-03 

3 . 91793E-03 

3. 92281E-03 

3 . 91793E-03 

3. 91304E-03 

3.91793E-03 

3.91304E-03 

3.90816E-03 

3.91304E-03 

3.90816E-03 

3.91793E-03 

3.91793E-03 

3.91304E-03 

3. 90816E-03 

3. 91793E-03 

3.91793E-03 

3.91304E-03 

3.91304E-03 

3.91304E-03 

3.90816E-03 


TIME 

« ft 




0 .2 .4 

,6 .8 



1 1.2 1.4 

1.6 1.8 




2 2.2 2.4 

2.6 2.8 




3 3.2 3.4 

3.6 3.8 




4 4.2 4.4 

4.6 4.8 




5 5.2 5.4 

5.6 5.8 




6 6.2 6.4 

6.6 6.8 




7 7.2 7.4 

7.6 7.8 




8 8.2 8.4 

8.6 8.6 




9 9.2 9.4 

9.6 9.8 





10 

11 

12 

10.2 

11.2 

12.2 

10.4 

11.4 

12.4 

10.6 

11.6 

12.6 

10.8 

11.8 

12.8 


13 

13.2 

13.4 

13.6 

13.8 


14 

14.2 

14.4 

14.6 

14.8 


15 

15.2 

15.4 

15.6 

15.8 


16 

16.2 

16.4 

16.6 

16.8 


17 

17.2 

17.4 

17.6 

17.8 


18 

18.2 

18.4 

18.6 

18.8 


19 

19.2 

19.4 

19.6 

19.8 


20 

20.2 

20.4 

20.6 

20.8 


~ir 

•2T75" 


”21' .3 


22 

22.2 

22.4 

22.6 

22.8 


23 

23.2 

23.4 

23.6 

23.8 

, 

24 

24.2 

24.4 

24.6 

24.8 


25 

25.2 

2574 1576 

""25.8 


26 

26.2 

26.4 

26.6 

26.8 


27 

27.2 

27.4 

27.6 

27.8 


28 

28.2 

28.4 

28.6 

28.8 


— zr 

~ir:r 

-zr.Tr 

27TF^ 2775 
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30 30.2 30.4 30.6 30. 8 

. li JU -JU* _31^ _3i A 

32.0001 32.2001 32.4001 32.6001 32.8001 

33.0001 33.2001 33.4001 33.6001 33.8001 

34.0001 34.2001 34.4001 34.6001 34.8001 


. UVU l 

36.0001 

i. 

36.2001 

36.4001 

36.6001 

36.8001 

37.0001 

37.2001 

37.4001 

37.6001 

37.8001 

38.0001 

38.2001 

38.4001 

38.6001 

38.8001 

39^00.1 _39.,.40fii _59^.09.1 . 

.3.9^0.01 

40.0001 

40.2001 

40.4001 

40.6001 

40.6001 

41.0001 

41.2001 

41.4001 

41.6001 

41.8001 

42.0001 

42.2001 

42.4001 

42.6001 

42.6001 

43.0001 

43.2001 

43.4001 

43.6001 . 

43.6001 

44.0001 

44.2001 

44.4001 

44.6001 

44.8001 

45.0001 

45.2001 

45.4001 

45.6001 

45.8001 

46.0001 

46.2001 

46.4001 

46.6001 

46.8001 

47.0001 

47.2001 

47.4001 

47.6001 

47.6001 

48.0001 

46.2001 

48.4001 

46.6001 


CYCLE * 5 

BEGTIME = 

599.6 

STRESS l * • 




15722.9 

15695.6 

15722.9 

15722.9 

15895.6 

15722.9 

15895.6 

15895.6 

15895.6 

15895.6 

15695.6 

15722.9 

15695.6 

15895.6 

15722.9 

15895.6 

15722.9 

15695.6 

15895.6 

15695.6 

15895.6 

15722.9 

15895.6 

15895.6 

15895.6 

16068.4 

15895.6 

15895.6 

15722.9 

15895.6 

15895.6 

15722.9 

15895.6 

15722.9 

15895.6 

15895.6 

15895.6 

15895.6 

15895.6 

15722.9 

15895.6 

15895.6 

15895.6 

15895.6 

15895.6 

15895.6 

15895.6 

15895.6 

15895.6 

15695.6 

15722.9 

15895.6 

15895.6 

15695.6 

15895.6 

15895.6 

16068.4 

15695.6 

15895.6 

15695.6 

15895.6 

15895.6 

15722.9 

15895.6 

16068.4 

15895.6 

16066.4 

15695.6 

15895.6 

15695.6 

15695.6 
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